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PREFACE 


This  report  is  the  second  of  a  sequence  of  two  reports  dealing  with  the  study 
of  possible  inertial  navigation  systems  for  submarine  use.  The  first  report, 
entitled  "Theoretical  Background  of  Inertial  Navigation  for  Submarines",  was 
issued  in  March  1951  and  dealt  with  the  fundamental  problem  of  inertial  navi¬ 
gation  with  the  discussion  restricted  to  idealized  systems.  The  present  report, 
the  second  of  this  sequence,  entitled  "Characteristics  of  Systems  Feasible  for 
Inertial  Navigation  of  Submarines"  considers  characteristics  of  nonideal  equip 
ment  that  must  be  used  in  any  practical  system.  On  this  basis,  one  system  out  oi 
the  many  theoretically  possible  types  has  been  chosen  as  most  nearly  meeting 
the  requirements  for  submarine  operation. 

This  work  has  been  carried  out  under  Contract  N50ri-07850  for  the  Office 
of  Naval  Research.  The  principal  responsibility  and  supervision  of  the  report 
has  been  carried  out  by  Mr.  Forrest  E.  Houston,  Assistant  Director  of  the 
Instrumentation  Laboratory,  with  the  assistance  of  Mr.  John  Hovorka.  Acknowl¬ 
edgment  is  also  due  the  Technical  Publications  Division  of  Jackson  A  Moreland 
for  the  preparation  of  the  report. 


Walter  Wrigley 

Cambridge,  Massachusetts 
October,  1951 
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ABSTRACT 


/  *  This  report  is  the  result  of  a  study  of  possible  submarine  inertial  navigation 

methods.  Prior  investigation  has  narrowed  the  field  of  inquiry  to  gravity-field 
navigation,  and  pointed  to  a  system  design  based  on  a  stable  vertical  element 
which  is  essentially  an  equivalent  pendulum  with  eighty-four-minute-period 
characteristics,  a  so-called  Schuler-tuned  system.  The  equivalent  pendulum  is 
represented  by  a  controlled  member  which  Is  oriented  by  accelerometer-monitored 
gyro  units  to  indicate  the  vertical.  The  accelerometer  units  furnish  data,  which, 
after  processing,  are  used  to  precess  the  gyro  units  and  drive  the  controlled 
member  so  that  it  responds  to  changes  in  the  vertical  direction  as  the  ship  moves 
along  its  course  over  the  surface  of  the  Earth.  The  gyro  units  also,  in  conjunc¬ 
tion  with  the  same  controlled  member  drives,  simultaneously  isolate  motions  of 
the  base  from  the  controlled  member.  This  kind  of  vertical  indicator  is  described 
in  detail  in  the  report,  with  the  discussion  including  such  problems  as  the  effects 
of  damping,  uncertainties  in  component  outputs,  and  inaccuracies  in  inputs  from 
external  sources.  With  this  system,  Schuler  tuning  eliminates  the  requirement 
for  precise  external  compensation.  The  optimization  of  system  parameters  is 
discussed,  and  numerical  values  are  given  for  optimum  parameters. 

When  the  indicated  vertical  has  been  established  by  a  controlled  member,  the 
problem  of  position  indication  can  be  solved  by  two  general  classes  of  systems. 

In  one  class,  accelerometer  data,  before  or  after  processing,  are  integratefckln  an 
open-chain  configuration  to  indicate  position  change  from  a  known  point  of  origin. 
The  accelerometer  itself  may  be  doubly-integrating  in  character,  in  which  case 
position  change  is  indicated  by  utilizing  its  output  directly.  This  class  of  position 
indicators,  deriving  position  information  from  integrators  appended  to  the  verti¬ 
cal-indication  loops  or  by  direct  double  integration  of  acceleration,  is  characterized 
by  simplicity  of  concept  and  geometry,  i.e.,  only  three  concentric  gimbals  are  re¬ 
quired,  no  more  than  are  required  for  vertical  indication  alone.  However,  the 
integrators,  being  open  chain  in  nature,  also  integrate  false  ground-speed  compo¬ 
nents  associated  particularly  with  gyro-unit  drift  and  w'th  Inaccuracy  in  Earth- 
rate  compensation.  This  leads  to  errors  in  indicated  position  which  are  cumulative, 
depending  on,  among  other  things,  the  time  of  operation  of  the  system. 

In  view  of  thp  long  operating  intervals  probably  required  of  a  submarine  navi¬ 
gation  system,  methods  for  avoiding  such  cumulative  errors  are  important,  in 
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this  report,  it  is  shown  that  position-indication  systems  of  this  class,  including 
those  using  a  doubly-integrating  accelerometer,  have  intrinsic  cumulative  errors 
in  their  output,  which  are  reducible  only  by  instrumentation  sufficiently  precise 
to  keep  ihe.se  errors  within  preassigned  limits  during  the  long  operating  intervals. 

The  second  class  of  position  indicators  discussed  requires  the  use  of  two 
additional  gimbals,  a  disadvantage  that  is  at  once  offset  by  the  absence  of  cumu¬ 
lative  error  in  latitude  indication,  and  by  efficient  minimization  of  cumulative 
error  in  longitude,  using  a  longitude  indicator  which  operates  in  conjunction  with 
the  latitude  indicator.  (Longitude,  because  of  its  arbitrary  nature,  lias  an  inherent 
possibility  lor  cumulative  error  in  its  indication:  this  is  not  true1  of  latitude.) 

This  class  of  systems  consists  of  two  general  types:  the  prealigned  system  and 
the  self-erecting  system.  Prealignment  of  ine  rtial  gyro  unit  axes,  so  that  they 
serve  as  "star-lines,"  presents  some  problems  in  a  submarine  installation,  par¬ 
ticularly  in  the  requirement  for  an  azimuth  reference.  A  self-erecting  system 
avoids  these  difficulties. 

The  self-erecting  latitude  indicator  discussed  in  this  report  operates  in  con¬ 
junction  with  an  azimuth  indicator  (an  inherent  gyrocompass)  and  with  the  verti¬ 
cal  indicator,  as  a  closcd-ioop  system  which  causes  a  reference  direction  on  a 
fourth  gimbal,  mounted  concentrically  with  the  controlled  member,  to  track  the 
Earth  polar  axis.  A  fifth,  and  innermost  gimbal,  mounted  within  the  latitude  gim- 
bal,  carries  an  integrating  gyro  unit  which  serves,  with  the  innermost  gimbal 
drive,  to  integrate  celestial  longitude  rate.  Indicated  longitudi  is  derived  from 
this  integrated  rate,  while  latitude  is  indicated  by  the  orientation  of  the  polar- 
axis-tracking  gimoal  relative  to  the  controlled  member.  This  position-indication 
method  is  examined  from  the  point  of  view  of  practical  instrumentation  and  accu¬ 
racy  of  indication.  The  problem  of  ground-speed  indication  with  this  system  is 
also  discussed,  and  it  is  shown  that  ground  speed  can  be  obtained  with  few  addi¬ 
tional  components. 

A  balancing  of  the  advantages  and  disadvantages  of  these  systems,  while  keep¬ 
ing  in  view  the  practicalities  of  instrument  design  and  the  demands  of  shipboard 
operation,  permits  the  recommendation  of  the  last- mentioned  position  indicator 
as  the  most  promising.  The  study  of  this  position  indicator  has  been  sufficiently 
exhaustive  to  permit  the  signment  of  numerical  values  to  the  important  system 
parameters.  These  are  ^ven  in  the  report.  In  addition,  potentially  profitable 
further  research  concerning  certain  specified  aspects  of  the  problem  is  suggested. 
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INTRODUCTION 


Inertial  guidance  nf  submarines  is  a  method  for  deriving  essential  navigation 
information  -  latitude,  longitude,  heading  and  speed  —  from  a  system  that  obtain: 
its  data  entirely  within  the  submarine.  Once  initial  conditions  have  been  inserted 
into  the  system,  it  must  operate  as  nearly  independently  of  external  information 
as  possible,  indicating  the  above  quantities  with  sufficient  precision  and  continuity 
for  the  successful  navigation  of  the  submarine. 

The  preceding  report  in  this  group,  Report  R-9,  Part  contains  a  discus¬ 

sion  of  fundamental  approaches  to  this  problem.  The  basic  method  discussed 
there  involves  a  stabilized  platform  as  an  essential  feature.  A  stabilization  servo 
system  receives  specific  force  data  from  accelerometers  or  pendulums,  and  de¬ 
livers  these  data  to  gyro  units  and  platform  drives,  which  cause  the  platform  to 
indicate  the  vertical.  Two  basic  methods  of  utilizing  the  indicated  vertical  are 
then  available  to  obtain  the  submarine's  essentially  geocentric  angle  of  travel. 

One  method,  referred  to  as  analytical  integratio  .  m  Part  I,  involves  the  integra¬ 
tion  of  accelerometer  outputs  or  of  gyro  inputs,  with  the  accelerometers  and  gyro 
units  mounted  on  the  stable  platform.  The  other  method,  called  geometric  inte¬ 
gration,  involves  the  direct  comparison  of  the  indicated  vertical  with  a  reference 
vertical.  The  conclusion  arrived  at  in  Part  1  was  that,  within  the  submarine,  a 
servo-driven  platform  on  which  were  mounted  two  accelerometers  (or  pendulums) 
and  two  gyroscopes  could  be  used  to  indicate  the  vertical,  and.  with  the  addition  of 
a  third  gyroscope,  to  indicate  heading.  The  vertical  indicator  serves  to  show  the 
direction  of  the  Earth's  gravity  field,  and  the  azimuth  indicator  shows  the  projec¬ 
tion  of  the  Earth's  polar  axis  in  the  horizontal  plane.  This  three -axis  stabilized 
platform  is  the  basis  for  various  methods  of  position  indication  described  in  the 
preceding  report,  and  in  addition  inherently  provides  data  on  heading,  roll,  pitch 
and  speed. 

It  is  the  function  of  this  report,  first,  to  examine  the  practicalities  of  instru¬ 
mentation  of  ?  platform  stabilized  about  the  vertical  and  in  azimuth,  and  second, 
to  re-examine  the  position-indication  methods  in  the  light  of  necessarily  imper¬ 
fect  instrumentation  and  imperfect  external  compensations.  This  discussion 
serves  for  the  rejection  of  systems  which,  while  they  appeared  theoretically 

*  Report  R  -9,  Part  I,  will  hereinafter  be  referred  to  as  Part  I.  Superscript  numbers  in  parentheses  lelr.i  lo  the  bibhogtaphy. 
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feasible  as  described  in  Part  I,  present  engineering  difficulties  in  execution. 

It  is  thus  possible  to  apply  certain  principles  oi  selection  and  compromise, 
which  are  discussed  in  the  following,  to  arrive  at  what  appears  to  be  the  most 
practical  submarine  ineriial  navigation  system.  The  discussion  of  this  system 
makes  up  the  largest  part  of  the  report.  Appended  to  the  report  is  a  series  of 
derivations  which  render  the  textual  discussion  in  mathematical  terms. 
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GENERAL  CONSIDERATIONS  CONCERNING  SUBMARINE 
INERTIAL  NAVIGATION 


The  essential  problem  to  be  discussed  is  the  determination  of  the  present 
position  of  the  submarine  on  the  Earth,  from  data  obtained  within  the  submarine. 
The  position  is  defined  in  terms  of  astronomical  latitude  and  longitude,  and  these 
are  in  turn  defined  as  two  geocentric  angles.  Astronomical  latitude  is  defined  as 
the  angle,  in  the  meridian  plane,  between  the  local  vertical  and  the  vertical  at  the 
equator.  Astronomical  longitude  is  defined  as  the  angle  between  the  meridian 
plane  containing  the  local  vertical  and  a  reference  meridian  plane.  It  is  to  be 
noted  that  astronomical  position  and  map  data  are  not  necessarily  identical.* 
Position  determination  thus  reduces  to 

1)  determination  of  the  local  vertical 

2)  determination  of  the  angle  between  the  local  vertical  and  two  Earth  refer¬ 
ence  directions. 

The  two  Earth  reference  directions  mentioned  above  may  be,  respectively,  a 
line  parallel  to  the  Earth  polar  axis  as  indicated  by  the  navigation  system,  and  a 
longitude  reference  direction  set  into  the  system  (latitude  and  longitude  are  then 
derived  directly  from  angles  between  the  indicated  directions);  or,  alternatively, 
the  original  vertical  and  azimuth  at  the  start  of  the  navigation  problem  may  serve 
as  references  (as  in  cases  where  rate  or  acceleration  integration  methods  are 
used  to  give  the  geocentric  angle  of  travel  of  the  submarine  from  the  starting 
point). 

In  Part  I,  ideal  components  formed  the  navigation  systems  considered.  These 
components  could  of  themselves  introduce  no  errors  into  a  system.  Practical  en¬ 
gineering  will  require  the  consideration  of  nonideal  behavior  of  components.  This 
nonideal  behavior  can  be  characterized  in  various  ways:  for  example,  the  over-all 
sensitivity  of  the  component  can  be  regarded  as  uncertain,  i.e.,  there  is  an  uncer¬ 
tainty  in  the  relation  of  the  output  to  the  input.  Or,  as  will  be  assumed  in  this  re¬ 
port,  the  component  sensitivities  may  be  regarded  as  pre-set  with  negligible  error, 
although  some  uncertainty  may  nevertheless  occur  in  the  output  of  a  component. 

*  Astronomical  position  data  are  derived  from  a  knowledge  of  the  local  vertical,  while  map  data,  which  are  in  effect  smoothed  astro¬ 
nomical  data,  ate  oetived  from  geodetic  triangulation  techniques.^'^'1^  The  corrections  required  to  bung  the  two  sets  of  position  data 
into  coincident  are  in  general  quite  small;  but  if  desired,  they  may  be  applied  in  a  manner  similar  to  that  used  to  make  variation  col¬ 
lections  to  a  magnetic  compass  reading. 
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In  any  case,  the  nonideal  character  of  a  component  will  he  displayed  ultimately  (in 
the  performance  equation)  by  an  uncertainty  term  appropriate  to  the  nature  of  the 
component  and  to  its  role  in  the  system.  The  word  uncertainty  here  refers  speci¬ 
fically  to  the  indicated  output  ot  the  component,  and  tin  indicated  quantity  uncer¬ 
tainty  is  defined  as  the  indicated  quantity  minus  some  *  tatistieal  average  of  this 
quantity.  Similarly,  the  indicated  quantity  error  is  defined  as  the  average  mdica- 
ed  quantity  minus  the  correct  value  of  the  quantity.  The  negative  of  the  error  is 
defined  as  the  correction  to  the  indicated  quantity.  I  inally,  the  indicated  quantity 
inaccuracy  is  defined  as  the  indicated  quantity  minus  the  correct  value  of  the  quan¬ 
tity.  Therefore,  at  any  instant,  the  inaccuracy  is  equal  to  the  error  plus  the  un¬ 
certainty  as  here  defined.  These  distinctions  are  profitable  because,  in  the  case 
of  a  component  within  the  system,  during  the  "transient"  stage  of  the  motion  the 
component  output  error  due  to  the  transient  regime  is  ultimately  correctible  by 
the  system  itself,  and  only  the  uncertainty  in  the  component  output  is  important. 
But,  in  the  case  of  quantities  brought  into  the  system  in  open-chain  configurations 
from  external  sources,  the  system  can  influence  neither  the  error  nor  the  uncer¬ 
tainty,  and  the  sum  of  these,  the  inaccuracy,  becomes  the  important  determinant 
of  system  behavior. 

DETERMINATION  OF  THE  LOCAL  VERTICAL 

The  definition  of  astronomical  position  given  above  requires  that  the  vertical 
at  the  present  position  of  the  submarine  be  indicated  by  the  navigation  system.  On 
a  stationary  base  the  vertical  is,  by  definition,  the  direction  along  which  a  plumb 
line  settles.  On  an  accelerating  base,  a  plumb  line  no  longer  indicates  the  true 
vertical;  instead,  it  tends  to  indicate  the  direction  of  the  apparent  vertical,  which 
is  a  function  of  the  acceleration  of  the  base.  It  was  pointed  out  in  Part  I  that  a 
useful  vertical  indicator  on  a  moving  base  can  be  conceived  of  as  an  accelero¬ 
meter-monitored  servo-controlled  platform  which  acts  as  a  lour  period  equiva¬ 
lent  pendulum.  The  dynamics  of  this  system  w’ould  be  adjusted  by  Schuler  tuning. 
This  tuning  consists  in  selecting  an  undamped  natural  resonant  frequency  for  the 
equivalent  pendulum,  consisting  of  the  platform  with  its  control  systems,  such 
that  the  error  in  *he  indication  of  the  vertical  by  the  platform  is  substantially  in¬ 
dependent  of  horizontal  accelerations  of  the  base. 

There  is  a  subtle,  basic  difference  between  this  kind  of  vertical  indicator  and 
conventional  stable  verticals  and  artificial  horizons.  The  latter  two  navigational 
aids  represent  a  philosophy  in  which  gravity  and  accelerations  are  considered  as 
essentially  different  physical  factors.  They  require  precise  compensation  inputs, 
that  is,  ship-motion  data  must  be  available  independently  and  of  the  same  order  of 
precision  as  that  desired  in  the  indication  of  the  vertical  direction.  This  is  due  to 
the  fact  that  the  vertical  is  associated  witli  gravity  alone,  and  accelerations  are 
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considered  as  interferences.  However,  the  Schuler-tuned  vertical  indicator  dis¬ 
cussed  above  does  not  require  precise  compensation  inputs  to  assure  small  steady - 
state  errors.  Fundamentally,  this  is  because  in  a  Schuler-tuned  system  the  indi¬ 
cated  vertical  is  not  associated  with  gravity  alone,  nor  are  accelerations  of  the 
base  to  be  regarded  as  interferences;  instead,  the  vertical  and  accelerations  of 
the  base  are  interrelated  through  the  geometry  and  kinematics  of  motion  over  the 
essentially  spherical  earth,  with  the  resultant  minimization  of  the  importance  of 
external  compensation  data.  It  is  true  that  compensation  inputs  w  ill  be  required 
by  a  practical  vertical  indicator  of  the  Schuler-tuned  type,  but  their  required  pre¬ 
cision  Is  such  that  complete  failure  of  compensations  should  not  seriously  impair 
the  system  performance;  they  act  as  a  factor  of  safety.  On  the  other  hand,  sys¬ 
tems  without  Schuler  tuning  cannot  achieve  required  performance  unless  the  com¬ 
pensation  is  of  a  high  order  of  precision. 

The  essential  operating  components  of  the  vertical  indicator  under  discussion 
are  shown  in  Fig.  1.  Tw’o  single-axis  accelerometers  (or  pendulum  units)  and 
three  single-axis*  integrating  gyro  units  are  shown.  This  system  is  described  on 
p.  22  of  Part  I.  In  terms  of  the  component  functions,  the  system  is  described  there 
from  the  standpoint  of  three  subsystems: 

1)  a  system  for  obtaining  data  on  the  resultant  specific  force  on  the  platform 
(controlled  member)  by  means  of  either  accelerometers  or  pendulums; 

2)  a  system  for  orienting  the  controlled  member  in  response  io  data  from  (1), 
so  that  the  controlled  member  indicates  the  vertical;  this  involves  gyro 
units  and  servo  drives; 

3)  a  system  for  modifying  data  from  (1)  before  they  are  applied  to  (2),  to  con¬ 
trol  the  dynamic  performance  of  the  whole  system. 

A  gyro  unit  is  characterized  by  an  ability  to  maintain  a  fixed  direction  in  iner¬ 
tial  space,  when  nj  torques  are  applied  tr  •«.  In  this  system,  torques  applied  by 
the  roll  and  pitch  of  the  ship  are  removed  oy  fast  gyro-monitored  servo  drives, 
which  orient  the  controlled  member  and  the  gyro  units  mounted  on  it  with  respect 
to  the  outer  gimbals;  and  at  the  same  time,  torques  applied  by  the  processed  ac¬ 
celerometer  data,  acting  to  process  the  gyro  units  more  slowly,  cause  the  control¬ 
led  member  to  rotate  approximately  geocentrically  at  the  angular  velocity  of  the 
indicated  vertical.  A  second  kind  of  system  is  conceivable,  in  which  the  gyro  units 
are  mounted  on  a  gimbal  separate  from  the  controlled  member  carrying  the  ac¬ 
celerometer  units,  in  which  case  the  controlled  member  can  be  oriented  with  re¬ 
spect  to  the  giinbal  carrying  the  gyro  units.  This  gimbal  then  remains  fixed  in  in¬ 
ertial  space,  serving  ultimately  as  a  source  of  reference  directions,  or  'star  lines", 

*  The  choice  of  single- axis  gyro  and  accelcromcfei  units,  as  conipaieil  with  two-axis  units,  has  hern  made  on  tup  hasis  ot  oarlio 
studies-^*. 
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Lin*  schematic  diagram  showing  essential  geometrical  parts  of  a  complete  verticol  indicating  system  using  two  single¬ 
axis  accelerometers  and  three  single-degree-of-freedom  gyro  units. 


carried  within  the  system,  for  position  indication.  Interference  isolation  against 
dpck  motion  is  thpn  obtained  hv  the  gyro-monitored  gimbnl  drives,  and  the  con 
trolled  member  drives  need  no  gyro  units.  The  practicability  of  this  scheme  in 
submarine  navigation  will  be  discussed  later  in  this  report.  In  both  this  system 
and  that  shown  in  Fig.  1,  the  principal  function  of  the  gyro  units  is  to  obtain  base 
motion  isolation,  and  the  angles  between  the  gimbals  external  to  the  controlled 
member  Indicate  roll  and  pitch  of  the  base,  and  therefore  of  the  submarine,  inci¬ 
dental  to  the  indication  of  the  vertical  by  the  controlled  member. 


DETERMINATION  OF  POSITION 

Once  the  problem  of  the  indication  of  the  local  vertical  has  been  solved,  an 
attempt  can  be  made  to  indicate  position  with  the  data  implicit  in  the  vertical  in¬ 
dicator.  For  example,  suppose  it  is  desired  to  indicate  latitude  by  integrating, 
with  respect  to  time,  the  input  signal  to  the  y-gyro  unit,i.e.,  the  unit  controlling 
motion  about  the  east-west  axis.  This  signal  is  proportional  to  the  y -component 
of  the  angular  velocity  of  the  indicated  vertical,  and  its  first  integral  is  approxi¬ 
mately  proportional  to  the  angular  displacement  of  the  indicated  vertical  in  the 
meridian  plane,  from  some  starting  position  which  serves  as  a  reference.  The  ap¬ 
proximation,  however,  depends  strongly  on  the  instrumentation,  and  can  lead  to 
serious  difficulties  tn  maintaining  correct  position  indication  at  all  points  on  the 
course. 

It  was  mentioned  before  that  two  classes  of  torques  act  to  precess  the  gyro 
units;  those  due  to  base  motion,  i.e.,  roll  and  pitch  of  the  ship,  and  those  due  to 
processed  accelerometer  data.  A  third  source  of  torques  on  a  gyro  unit  lies  in 
such  phenomena  as  mass -unbalance  and  friction  within  the  gyro  unit  itself.  These 
will  be  referred  to  as  uncertainty  torques.  They  give  rise  to  an  uncertainly  in  the 
gyro  unit  gimbal  angular  velocity  called  gyro  unit  drift.  This  drift  is  indistinguish¬ 
able  in  the  gyro  unit  output  from  the  true  angular  velocity  of  the  indicated  vertical; 
i.e.,  gyro  drift  produces  the  effect  of  a  false  ground  speed  component. 

There  is  a  similar  additional  source  of  false  ground  speed  which  comes  about 
in  the  following  way:  Suppose  that,  by  means  of  a  gyro  compass  operating  in  con¬ 
junction  with  the  z-gyro  in  Fig.  I,  the  controlled  member  is  oriented  with  the  x- 
axis  tracking  true  north,  so  that  the  y-axis  is  pointed  approximately  east.  Then 
the  x-axis  gyro  unit  will  sense  the  horizontal  projection  of  the  sidereal  rotation  of 
the  Earth  (hereinafter  referred  to  as  horizontal  Earth  rale),  while  the  y-axis  gyro 
unit  will  not.  The  x-axis  gyro  unit,  if  the  controlled  member  is  to  be  kept  in  the 
horizontal  plane,  must  be  compensated  for  horizontal  Earth  rate,  i.e.,  this  parti¬ 
cular  input  angular  velocity  component  must  be  cancelled  in  some  way.  One  pos- 
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sibility  is  to  prwoss  the  gyro  unit  with  a  signal  proportional  to  horizontal  Earth 
rate.  This  result;,  in  a  compensation  torque  on  the  gyro  unit  which  exists  in  ad¬ 
dition  to  the  torques  discussed  previously.  Lack  of  precision  in  this  Earth-rate 
compensation  input  will  result  in  a  residual  gyro  unit  gimbal  angular  velocity  in¬ 
accuracy.  This  inaccuracy  is  functionally  identical  with  the  aforementioned  gyro 
unit  drill  in  its  effect;  combined,  these  constitute  a  false  ground  speed  component. 

Note  that,  since  these  errors  are  essentially  dynamic,  they  leave  the  stead\ 
state  indicated  vertical  unattected.  This  i  so  because  the  system  is  accelerometer 
monitored,  and  the  integrator  between  the  accelerometer  and  the  gyro  automatically 
acquires  the  proper  bias  to  null  this  false  ground  speed  component;  the  indicated 
vertical  is  unaltered.  These  phenomena  are  common  to  all  systems  which  use  ac¬ 
celerometers  or  pendulums  to  track  the  apparent  vertical  and  w  hich  use  processed 
gyro  units  to  effect  base  motion  isolation. 

The  aforementioned  false  ground  speed  components  are  necessarily  included 
in  the  integration  of  rate  or  acceleration  data  to  indicate  position,  and  lead  to  an 
error  in  position  which  increases  with  time,  further  examination  will  show  that 
such  cumulative  errors  are  unavoidable  with  this  method  of  indicating  either  lati¬ 
tude  or  longitude. 

Similarly,  the  direct  double  integration  of  the  angular  acceleration  of  the  in¬ 
dicated  vertical,  by  whatever  method  -  e.g.,  a  doubly- integrating  accelerometer  - 
involves  again  a  cumulative  error  in  the  indicated  position.  A  false  acceleration 
component,  corresponding  to  the  lalse  ground  speed  component  in  another  part  of 
the  vertical  indicator  loop,  is  responsible  for  the  cumulative  error  in  this  case. 
Thus,  these  "open-chain  integration  and  direct  double -integration  methods  of 
indicating  position,  while  having  the  virtue  of  relative  simplicity,  simultaneously 
present  serious  drawbacks  in  practice. 

A  further  difficulty  enters  this  problem  when  the  possibility  of  integrator  drift, 
a  common  characteristic  of  such  devices,  is  considered.  In  the  accelerometer  data 
processing,  drift  in  the  integrators  used  will  affect  the  indicated  ground  speed, 
while  in  open-chain  integrators  used  to  determine  position,  integrator  drift  will 
lead  to  cumulative  errors  in  the  indicated  position.  The  eniire  problem  of 
cumulative  errors  is  traceable  to  the  fact  that  only  indicated,  not  true,  motions 
are  available  as  data  from  the  vertical  indicator;  more  generally,  aboard  the  ship, 
motions  of  the  various  parts  of  the  system  with  respect  to  each  other  and  to  the 
ship  are  the  only  motions  from  which  the  position  of  the  ship  can  be  deduced. 

It  is  therefore  of  interest  to  consider  the  possibility  of  establishing  an  inertial 
reference,  fixed  in  space,  and  relatively  independent  of  the  accelerometer-moni¬ 
tored  controlled  member.  As  was  mentioned  in  the  discussion  of  vertical  indica¬ 
tion,  the  accelerometers  of  f  ig.  1  might  be  mounted  on  the  controlled  member  as 
shown,  but  with  the  three  gyro  units  mounted  on  a  separate  platform,  which  re¬ 
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mains  fixed  in  inertial  space;  then  the  controlled  member  can  be  rotated  with  re¬ 
spect  to  the  platform  by  a  sidereal  time  drive,  about  a  polar  axis"  established  on 
the  gyro  untt  platform.  Evidently,  latitude  and  longitude  can,  with  suitable  initial 
gyro  unit  alignment,  be  indicated  directly  as  angles  between  gimbals.  Only  two 
cumulative  errors  are  then  possible:  that  due  to  gyro  unit  drift,  and  that  due  to 
inaccuracy  in  the  sidereal  time  drive.  At  the  present  time  accuracies  of  the 
order  of  one  part  in  one  hundred  thousand  appear  attainable  in  matching  si¬ 
dereal  time  drive  angulur  velocity  to  Earth  rate.  It  may,  therefore,  be  expected 
that  the  gyro  unit  drift  will  be  the  major  source  of  error  in  a  system  of  this 
type. 

This  situation  appears  to  be  unavoidable  in  the  case  of  longitude  indication. 
Astronomical  longitude  is  nut  a  unique  quantity;  only  longitude  difference  with  re¬ 
spect  to  a  present  reference  can  lie  indicated.  The  indicated  longitude  w  ill  in  no 
case  be  more  precise  than  this  reference  setting.  Furthermore,  it  must  be  ob¬ 
tained  by  some  form  of  integration  (as  in  the  case  just  described,  where  the  gyros 
and  drives  act  as  integrators),  with  the  concomitant  cumulative  error. 

Latitude  indication  is  not  necessarily  bound  by  these  restrictions.  A  self- 
erecting  latitude  system  is  conceivable,  with  a  gimbal  which  moves  with  respect 
to  the  controlled  member  until  the  gimbal  is  parallel  to  the  Earth  polar  axis. 
Thereafter  the  gimbal  moves  at  the  indicated  latitude  rate,  about  the  east-west 
axis,  with  respect  to  the  controlled  member;  it  is  monitored  by  the  azimuth  stabi¬ 
lizer,  and  ultimately  by  the  y-axis  accelerometer.  The  angle  between  this  gimbal 
and  the  indicated  horizontal  is  then  the  indicated  latitude. 

These  matters  will  each  be  discussed  in  more  detail  on  the  pages  follow  ing. 


THE  DAMPED  VERTICAL  INDICATOR 


An  undamped  vertical  indicator  of  the  type  considered  in  Part  1,  pp.  28  if.,  is 
necessarily  subject  to  a  fixed-amplitude  continuous  oscillation  with  an  84 -minute 
period.  Such  a  system  might  be  practical  under  the  following  conditions: 

1)  The  time  of  operation  is  sufficiently  short,  i.e.,  of  the  order  of  one  period. 

2)  The  initial  conditions  are  set  into  the  system  with  sufficient  precision,  i.e,, 
such  that  the  amplitude  of  oscillation  plus  the  steady-state  error,  where  it 
occurs,  will  not  exceed  the  preass igued  error  in  the  indication  of  the  verti¬ 
cal. 

In  submarine  operation,  an  operation  time  of  the  order  of  at  least  eight  Schuler 
periods  must  be  allowed  for,  and  in  addition  the  precision  of  the  initial  settings  is 
likely,  for  practical  reasons,  to  be  somewhat  poorer  than  the  ultimate  accuracy 
desired  In  vertical  indication.  The  long  operating  time  allows  a  large  probability 
for  the  system  to  be  disturbed  during  a  run,  making  it  possible  to  incur  ei  rors 
larger  than  those  associated  with  uncaging  of  the  instrument.  Insufficient  preci¬ 
sion  in  initial -condition  settings  requires  that  the  system  reduce  its  amplitude  of 
oscillation  after  the  submarine  is  set  in  motion. 

Evidently,  an  undamped  system  satisfies  neither  requirement.  On  the  other 
hand,  a  damped  system  suffers  from  two  serious  limitations: 

1)  Insofar  as  damping  suppresses  high  input  frequencies,  it  lengthens  the 
solution  time.*  This  effect  is  important  both  when  the  system  is  initially 
uncaged,  and  when  it  suffers  transient  disturbances  while  the  submarine 
is  in  motion. 

2)  The  delay  in  achieving  a  solution  is  operative  at  other  times  also;  effec¬ 
tively,  a  damped  system  responds  to  values  of  the  input  averaged  over  a 
time  of  the  order  of  the  solution  time,  rather  than  to  instantaneous  input 
values.  The  result  is  manifested  as  errors  which  depend  on  the  rapidity 
of  the  inputs.  These  are  called  forced  dynamic  errors  (cf  Part  1,  p.  17). 


*  Solution  lime  is  rielmed  for  these  reports  as  !iie  time  reiiinrert  lor  a  system  to  achifpp  ninptv-'iv*>  percent  o'  the  change  a1"  ateit 
with  a  given  input  change. 
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This  unavoidable  compromise  between  short  solution  time  and  small  forced  dyna¬ 
mic  errors  confronts  the  designer  of  any  practical  system. 

Consider  the  vertical  indicator  of  Fig.  1.  The  simplest  damping  method  might 
involve  putting  a  direct  by -pass  around  the  integrator,  between  the  accelerometer 
and  the  gyro  unit.  However,  this  corresponds  to  a  high-pass  filter ,  ind  would  re¬ 
sult  in  a  system  that  is  excessively  responsive  to  high*  frequencies,  some  of  which, 
in  view  of  the  practicalities  of  instrumentation,  can  be  expected  to  be  spurious. 

The  resultant  system  is  analogous  to  a  pendulum  with  frictional  damping,  and  is 
describable  by  a  second-order  differential  equation. 

High-frequency  suppression  can  bo  obtained  by  combining  the  integrator  and 
by-pass  with  a  low -pass  filter,  which  can  be  represented  as  an  integrator  with 
direct  feedback.  This  kind  of  system  is  shown  in  Fig.  2.  Its  performance  equation 
is  derived  in  Derivation  1A.  (Derivations  are  included  in  the  Appendix,  beginning 
on  page  52).  The  differential  performance  equation,  Eq.  (1-16),  is  of  third  order  in 
time  derivatives  of  the  correction  to  the  indicated  vertical.  The  tuning  of  the  sys¬ 
tem  to  an  approximate  84-minute  period  (Schuler  tuning)  is  therefore  possible  in 
alternative  ways.  In  terms  of  the  physical  inputs  to  the  system,  Schuler  tuning 
may  be  accomplished  by: 

1)  minimizing  the  acceleration  error,  or 

2)  minimizing  the  jerk  (acceleration  derivative)  error. 

The  result  of  each  method  is  displayed  in  the  two  amalitude- ratio  vs  fre¬ 
quency-ratio  logarithmic  plots.  Fig.  3.  It  should  be  noted  that,  were  the  damping 
sufficiently  reduced  in  either  system,  the  system  would  approach  an  equivalent 
Schuler  pendulum  with  a  natural  frequency  corresponding  to  a  frequency  ratio  of 
unity  (since  the  reference  frequency  is  the  Schuler  frequency).  But  the  frequency- 
ratio  for  maximum  amplitude  ratio  is  higher  for  acceleration-error  minimization 
and  lower  for  jerk-error  minimization,  the  departure  of  the  frequency  ratio  from 
unity  being  determined  in  both  cases  by  the  damping.**  The  relative  heights  of 
the  peaks  in  the  curve  are  also  determined  by  the  damping. 

The  choice  between  the  two  response  functions  can  then  be  based  on  whether  it 
is  desirable,  relatively,  to  suppress  low  frequencies  or  high  frequencies;  i.e., 
whether  the  expected  forced  errors  due  to  necessarily  Imperfect  compensation  in¬ 
puts  will  be  primarily  of  a  high-  or  low-frequency  character.  (Instrumentation  er¬ 
ror  will  presumably  be  sufficiently  reducible  to  give  the  desired  over -all  accuracy). 
The  assumption  is  here  made  that  these  uncontrolled  forced  errors  will  be  short- 
period,  compared  with  84  minutes,  and  in  Derivation  1A,  the  jerk  error  is  minimized 


•  It  is  assumed  that  the  sjstem  in  this  discussion  will  ee  SohuteMuned.  so  that  the  letetemc  *'»auentv  corresponds  to  a  period  ol 
M  minutes. 

■*  Specifically,  by  the  dumping  ratio,  which  is  defined  a.,  the  ratio  of  the  actuai  damping  tn  rutical  damping. 
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Fig,  3.  Amplitude  ratio  (correction  to  indicated  verticol/ii.occurocy  in  water-speed 
compensation)  vs.  frequency  rotio  (input  frequeney/Schuler  frequency)  far  two  vertical 
indicators:  (1)  with  jerk  error  minimized  by  Schuler  tuning, acceleration  error  minimized 
by  choice  of  undomped  naturol  frequency,  ond  (2)  with  acceleration  error  minimized  by 
Schuler  tuning,  jerk  error  minimized  by  choice  of  undamped  notural  frequency.  Damping 

ratio  is  0.3  in  both  coses. 

by  an  adjustment  of  the  system  parameters, which  corresponds  to  Schuler  tuning 
(Eq.  (1-13)). 

The  acceleration  error,  however,  remains.  It  is  proposed  to  minimize  this 
with  a  signal  representing  the  appropriate  ship  velocity  component.  In  the  case  of 
the  x-system,  this  component  is  the  longitude  rate;  in  the  y -system,  the  latitude 
rate.  This  velocity  component,  hereinafter  referred  to  as  the  water-speed  com¬ 
pensation,  is  fed  from  die  ship  pitomeler  log,  through  the  azimuth  resolver,  to  the 
high-frequency  integrator;  the  effect  of  this  treatment  of  the  ship  velocity  signal 
is  to  minimize  the  acceleration  error,  if  the  signal  is  suitably  adjusted  (Eqs.  (1-14) 
and  (1-15)).  Actually,  the  error  could  be  made  more  nearly  zero  it  the  ship  velo¬ 
city  with  respect  to  the  Earth  were  available  as  a  datum;  even  a  perfect  pitometer 
log  can  give  only  the  water  speed.  Therefore,  ocean  current  gradients  (although 
not  the  ocean  current  itself)  are  included  along  with  uncertainties  in  the  water- 
speed  indication  in  the  velocity-compensation  inaccuracy  (see  Eq.  (1-14)). 
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The  differential  performance  equation,  Eq,  (1-16),  is, as  observed  third  order, 
with  the  following  principal  forcing  terms,  shown  on  the  right-hand  side  of  the 
equation  in  the  order  given: 

1)  A  term  proportional  to  the  rate  of  change  of  the  inaccuracy  of  the  water  - 
speed  compensation  with  respect  to  true  ground  speed. 

2)  A  term  proportional  to  the  rate  of  change,  and  second  rate  of  change,  of 
the  inaccuracy  in  the  Earth-rate  compensation  current,  and  a  gyro  drift 
error  term  proportional  to  the  rate  of  change,  and  second  rate  of  change, 
of  the  gyro  gimbal  angular  velocity  uncertainty. 

3)  A  term  proportional  to  the  accelerometer  output  uncertainty  and  to  its 
rate  of  change,  and  a  term  proportional  to  the  difference  between  the  rate 
of  change  of  the  uncertainty  in  the  output  of  the  low-frequency  integrator 
by-pass  and  the  rate  of  change  of  the  uncertainty  in  the  feedback  channel 
output.  These  constitute  an  acceleration  error. 

The  left-hand  side  of  £q.  (1-16)  is  concerned  with  the  system  dynamic  char¬ 
acteristics.  Derivation  IB  deals  with  the  assignment  of  numerical  values  to  these 
characteristics.  The  dynamic  performance  can  be  described  in  terms  of  the  damp¬ 
ing  ratio  and  the  undamped  natural  frequency  of  the  quadratic  term,  and  the  char¬ 
acteristic  time  of  the  first-order  term.  The  numerical  values  of  these  three  quan¬ 
tities  ultimately  determine  the  actual  values  of  component  sensitivities  (Eq.  (1-17)). 
These  three  parameters  are  chosen  for  convenience  in  the  system  design,  since 
they  are  universally  applicable  to  linear  systems  of  this  general  nature.  In  Deri¬ 
vation  IB,  the  numerical  values  are  assigned  according  to  the  following  procedures: 

1)  The  damping  ratio,  as  has  been  observed,  has  an  optimum  value  of  zero, 
but  the  resulting  undamped  system  would  be  objectionable  for  the  reasons 
given.  The  chosen  damping  ratio  is  therefore  a  compromise.  The  smallest 
damping  ratio  consistent  with  a  reasonable  solution  time  is  considered  to 
be  about  0.3. 

2)  The  resonant  frequency  for  the  undamped  system  is  selected  to  minimize 
the  acceleration  error  in  Eq.  (1-12),  before  the  water-speed  compensation 
is  applied.  It  is  desirable  to  do  this  because  of  the  aforementioned  ex¬ 
pected  difficulty  in  providing  a  perfect  velocity-input  in  making  the  com¬ 
pensation.  The  resonant  frequency  thus  found  is  determined  by  the  chosen 
damping  ratio,  and  is  95,4  percent  of  the  Schuler  frequency. 

3)  Given  the  damping  ratio  and  the  resonant  frequency  for  the  undamped  sys¬ 
tem,  the  characteristic  time  is  determined  by  Eq.  (1-24)  to  be  1.43  hours. 

The  transient  behavior  of  the  system  is  studied  in  Derivation  1C  by  consider¬ 
ing  the  homogeneous  equation,  Eq.  (1-27),  derived  from  Eq.  (1-16).  As  in  Deriva¬ 
tion  IB,  it  is  assumed  that  the  left-hand  side  of  Eq.  (1-16)  is  factorable  into  a  first- 
order  term,  containing  the  characteristic  time,  and  a  second-order  term,  containing 
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the  damping  ratio  and  the  natural  frequency.  Laplace  transform  methods  then 
furnish  the  solution  (Eq.  (1-30)).  The  correction  to  the  indicated  vertical  contains 
the  following  two  important  groups  of  terms: 

1)  A  damped  exponential,  controlled  by  the  characteristic  time. 

2)  A  damped  sinusoid,  with  a  maximum-amplitude  frequency  dependent  on  the 
damping,  as  has  been  observed;  the  frequency  ratio  at  the  frequency  of 
maximum  amplitude  is  equal  to  unity  minus  the  square  of  ihe  damping 
ratio.  The  damping  of  the  sinusoidal  term  is  controlled  by  the  damping  - 
ratio- naturrl-frequency  product. 

Refer  again  to  Eq.  (1-16)  of  Derivation  1.  On  the  right-hand  side  of  the  equa • 
tion  are  several  forced -error  terms,  which  w  ill  be  discussed  in  terms  of  a  method 
to  be  outlined  below.  The  most  important  forced  errors  are  caused  by  the  follow¬ 
ing: 

1)  velocity  compensation  inaccuracy, 

2)  gyro  unit  drift, 

3)  accelerometer  output  uncertainty. 

The  functional  form  of  these  errors  is  not  precisely  known;  this  obviates  the  pos¬ 
sibility  of  obtaining  a  particular  solution  to  Eq.  (1-16)  by  such  means  as  Laplace 
transform  methods.  However,  it  is  possible  to  exclude  certain  regions  from  the 
frequency  spectrum  of  each  of  these  forced  errors  as  being  highly  improbable. 
This  fact,  coupled  with  the  frequency  response  of  the  system  to  sinusoidal  varia¬ 
tions  in  each  of  the  forced  errors  taken  individually,  permits  some  statements  to 
be  made  concerning  the  effect  of  these  errors  on  the  over-all  accuracy  in  the  in¬ 
dication  of  the  vertical. 

Thus,  since  the  exact  frequency  spectrum  for  each  error  is  unknown,  it  can 
first  be  assumed  that  ajl_  frequencies  occur  with  equal  amplitudes,  i.e.,  that  the 
input  to  the  system  consists  of  "white  noise”.  A  harmonic  analysis  of  the  output 
then  gives  the  normalized  frequency  response.  Figure  4  shows  the  result  of  such 
an  analysis  applied  to  the  velocity-compensation  inaccuracy.  The  curve  exhibits 

a  peak  near  the  Schuler  frequency.  At  high  and  low  frequencies  the  response  is 
much  attenuated.  The  conclusion  is  that  velocity-compensation- inaccuracy  oscil¬ 
lations  having  frequency  components  near  the  Schuler  frequency  wilt  be  most  im¬ 
portant,  and  that  other  frequencies  will  be  less  so.  To  evaluate  the  net  effect  of 
the  velocity-compensation  inaccuracy  on  the  over-all  error,  then,  it  must  be 
known  whether  or  not  frequency  components  in  this  forced  error  are  near  the 
Schuler  frequency.  Little  data  are  available  on  this  subject,  but  the  presence  of 
84-minute  period  forcing  functions  for  an  appreciable  length  of  time  does  not  ap¬ 
pear  lik’ly. 
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Similar  considerations  apply  to  Kip-ares  5  and  G,  which  apply  to  the  gyro  drilt 
and  accelerometer  uncertainty  errors,  respectively. 

Since  this  is  a  single-axis  system,  and  two  such  systems,  geometrically  ortho¬ 
gonal,  are  required  for  vertical  indication,  the  coupling  between  two  similar  sys¬ 
tems  was  investigated,  as  follows:  An  equivalent  complete  vertical  indicator,  rep¬ 
resented  by  Eqs.  (1-39)  and  (1-40)  in  Derivation  ID,  was  set  up  on  a  Reeves 
Electronic  Analogue  Computer  (RF.AC)  in  the  Instrumentation  Laboratory 
The  simulated  vertical  indicator  was  effectively  uncaged  with  equal  errors 
in  the  vertical  about  both  the  x-  and  y-axes.  The  system  came  to  rest  in  the 
manner  of  a  damped  Foucault  pendulum -^(the  damping  ratio  chosen  was  0.3).  Pig- 
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Fig.  4.  Effect  of  water-speed  compensation  inaccuracy  on  vertical  indication,  os  a  function 

of  frequency. 
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uro  7b  shows  a  plot  of  x-  vs  y-  displacements  from  the  true  vertical  as  a  solution 
was  reached.  It  is  seen  that  the  projections  of  the  two-dimensional  oscillations  on 
the  x-  and  y-  input  planes  are  practically  the  same,  a  situation  that  would  occur 
for  decoupled  x-  and  y-  systems  (Kip.  7a).  The  conclusion  is  that  this  two-dimen¬ 
sional  pendulum  effect  can  be  neglected,  and  single-axis  results  are  immediately 
extensible  to  the  three-dimensional  case.  These  remarks  are  important  bec  ause 
in  subsequent  discussions  it  will  be  assumed  tliat  as  a  good  approximation,  the  x- 
and  y -systems  are  independent. 
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Fig.  5. Effect  of  gyro  unit  drift  on  verticol  indication  ,  os  o  function  of  frequency 
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Fig.  6-  Effect  of  occelerometer  uncertainty  on  vertical  indication,  as  a  function  of  frequency. 
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POSITION  INDICATION 


POSITION  INDICATION  BY  ANALYTICAL  INTEGRATION  CP  RATES 

A  distinction  whs  made  in  Part  1  (p.  27)  between  two  general  method:  of  poM 
tion  indieau  m:  analytic  integration  of  the  acceleronu  ter  output  data,  and  pcoii.e- 
tric  integration  of  the  accelerometer  output  data.  These  meth  ids  will  be  del  tiled 
again  here  in  the  same  way.  Analytu  1  integration  is  characterized  by  a  modifi¬ 
cation  of  the  accelerometer  output  data  by  means  of  two  integrators,  so  that  the 
Output  of  the  second  integrator  is  a  signal  proportional  to  the  essentially  geocen¬ 
tric  angular  displacement.  Geometric  integration,  on  the  othei  hand,  is  character¬ 
ized  by  a  physical  comparison  of  the  indicated  vertical  with  a  reference  vertical; 
the  angle  between  these  two  verticals  gives  the  indicated  position. 

First,  analytical-integrating  systems  will  be  discussed.  They  will  be  treated 
in  terms  of  single-axis  damped  systems,  with  vertical  indication  pr  ivtded  in  each 
case  by  a  system  similar  to  that  of  Fig.  2. 

Open-chain  integration  of  the  angular  velocity  of  the  indicated  vertical 

This  system  is  shown  in  Fig.  8.  A  performance  equation,  Eq,  (2-9).  is  derived 
in  Derivation  2A,  to  assess  the  effect  of  component  uncertainties  and  input  inac¬ 
curacies  (which  are  treated  as  in  Derivation  1).  The  correction  to  the  indicated 
position  as  given  by  Eq.  (2-9)  is  equal  to  the  following  sum: 

1)  the  initial  position  correction, 

2)  the  change  in  the  correction  to  the  indicated  vertical, 

3)  an  indicated  ground  speed  integral  with  the  time  of  operation  as  its  limits. 
The  most  serious  source  of  error  is  the  last  term,  in  which  F'arth-rate  cotttpi  n:-a - 
tion  inaccuracy,  gyro  drift,  and  the  position  integrator  drift,  w  hich  appear  as  false 
ground-speed  components,  are  integrated  with  respect  „o  time.  Tins  cumulative 
error,  unlimitt'd  except  by  the  time  of  operation,  is  the  direct  result  of  the  open- 
chain  nature  of  the  last  integration  giving  the  position.  Note  that  a  one -nautical - 
mile-per-hour  false  ground  speed  component  gives  rise  to  a  minute  of  arc  error 
in  the  indicated  position  for  each  hour  of  operation  time.  \  false  ground-speed 
component  of  this  magnitude  would  correspond  to  a  gyro  drift  of  about  0.015  de¬ 
gree  per  hour,  or  to  an  Earth-rate  compensation  signal  which  failed  to  be  constant 
to  within  less  than  0.1  percent  or  to  a  position-integrator  drift  of  om  minute  of 
arc  per  hour.  Several  hours’  operation  time  can  thus  be  expected  to  give  rise  to 
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excessive  errors,  unless  this  Ial*-e  (•.round  speed  compoiu  nt  car.  be  Kept  sufficiently 
small.  The  relative  simplicity  ot  conception  ol  th»s  method  ni  position  indication 
is  thus  offset  by  a  disadvantage  in  ilie  hi|'H  pr<  ci,  ion  of  the  instrumentation  rt 
quin  d  to  minimise  the  false  (‘round  spied  component.  The  magnitude  of  the  op¬ 
eration  time  is  a  strong  determining  lac  tor  in  judging  the  practicability  of  this 
system,  and  is  in  be  balanced  against  possible  instrumentation  problems,  which 
must  include  stability  anti  maintenance  con;  tderatiotis  m  the  initial  design  proc  e  ¬ 
dure.  In  submarine  navigation,  operation  time  intervals  ol  the  orde  r  of  several 
hours  are  required;  moreover,  the  operation  time  is  necessarily  somewhat  indefi¬ 
nite  in  any  case,  a  situation  to  be  contrasted,  for  example,  with  the  guidance  ol  a 
vehicle  over  a  preselected  course  between  a  given  starting  point  and  destination. 
Systems  of  this  type,  i.e.,  those  .subject  to  cumulative  errors,  are  then. lot  t  t  >  be 
regarded  with  disfavor  as  solutions  to  the  submarine  navigation  problem. 

Open-chain  double  integration  of  the  angular  acceleration  of  the  indicated  v«  meal 

This  system  is  shown  in  Tip.  9,  and  a  performance  equation.  Kq.  <2-37  t,  ts 
derived  in  Derivation  2B.  The  most  important  effect  of  component  uncertainties 
and  input  inaccuracies  is  to  create  two  cumulative  error  terms.  In  vi<  w  of  tin 
previous  discussion,  these  terms  will  now  be  examined.  One  of  the  terms  is  a 
single  integration,  with  respect  to  the  time  of  operation,  of  a  false  ground  sjieed 
component;  the  oilier  is  a  double  integration  of  a  false  acceleration  component. 

The  presence  of  a  term  of  this  form  is  traceable  to  the  open-chain  double  integra¬ 
tion  peculiar  to  this  method  of  position  indication.  In  connection  with  the  false 
ground  speed  component,  note  that  there  are  four  principal  sources  of  error: 

1}  a  term  proportional  to  the  initial  value  of  the  correction  to  the  indicated 
vertical, 

2)  gyro  drift, 

3)  Earth-rate  compensation  inaccuracy, 

4)  velocity  compensation  inaccuracy. 

A  false  ground  speed  component  of  one  nautical  mile  per  hour,  giving  rise  to  an 
error  of  one  minute  of  are  in  the  indicated  position  for  each  hour  nf  opera ti  n  time, 
corresponds  to  an  initial  correction  to  the  indicated  vertical  c.f  0.03  minute  ol  arc. 
or  to  a  gyro  drift  of  0.003  degree  per  hour.  Evidently,  the  difficulties  in  accurate 
position  indication  associated  with  cumulative  error  in  the  previous  ease  discussed 
are  here  even  more  serious.  Further  examination  of  the  error;  incurred  in  using 
an  open-chain  double-integration  method  confirms  this.  Therefore,  the  considera¬ 
tions  of  the  effect  of  cumulative  errors  on  a  submarine  navigation  system  given  in 
the  previous  rase  warrant  rejection  of  the  present  method  as  infeasible. 
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Direct  double  integration  of  .lccelerntti 'll 

This  system  is  shown  in  Tig.  10.  The  aeceleroineter  here  is  inherently  doubly 
integrating  (Part  I,  p.  34',  and  therelnre  is  an  iniinediate  source  ol  position  data. 

A  modification  of  ti>  •  vertical  indicator  ol  hip.  2  i:  used;  this  is  cl  is  cussed  at  the 
beginning  ol  Derivation  2C.  In  Derivation  2C  a  perlorinance  equation  (Kq.  ( 2 - G 7  i 
for  position  indication  is  also  derived.  The  most  significant  errors  are  cumulative 
as  in  the  two  previous  e  rnes,  i.i*„  proportional  to  a  la  1st*  ground  speed  component 
integrated  over  the  time  of  operation.  The  integrand  contains  the  lollowtng  two 
principal  sources  of  error: 

1)  gyro  drift, 

2)  Earth-rate  compensation  inaccuracy. 

The  same  conclusions  regarding  indicated  position  errors  apply  to  this  case  as; 
apply  to  the  case  of  the  single  open-chain  integration  method  already  discussed. 
That  is,  a  false  ground  speed  component  of  one  nautical  mile,  creating  a  minute  of 
arc  error  in  the  indicated  position  lor  each  hour  of  operating  time,  corresponds  to 
a  gyro  drift  of  about  0.015  degree  per  hour,  or  to  an  Earth-rate  compensation  er¬ 
ror  of  0.1  percent. 

Again,  although  the  method  appears  superficially  simple  and  therelore  of  po¬ 
tential  applicability,  the  long  time  of  operation  encountered  in  submarine  naviga¬ 
tion  would  probably  make  instrumentation  problems  difficult. 

Note  that  fhere  is  a  duplication  of  effort,  so  to  speak,  involved  in  differentiat¬ 
ing  the  already -integrated  position  data  to  obtain  the  velocity  signal  required  to 
precess  the  integrating  gyro  unit.  This  situation  can  be  avoided  by  using  a  rate- 
gyro  unit  and  omitting  the  differentiating  network,  but  this  would  involve  poorer 
drift  characteristics  in  the  gyro  unit  than  can  be  obtained  with  an  integrating  gyro. 
Retention  ot  the  differentiating  network,  on  the  other  hand,  involves  accentuated 
system  response  to  high  frequencies.  It  is  possible  that  this  problem  might  be 
resolved  by  further  study.  However,  the  cumulative  errors  would  still  be  present, 
and  appear  to  be  unavoidable.  For  this  reason  this  system  is  rejected  as  not  being 
applicable  to  the  submarine  problem. 

AZIMUTH  STABILIZATION 

All  of  the  foregoing  systems  require  some  form  of  azimuth  stabilization,  l.e.. 
orientation  of  the  controlled  member  in  the  horizontal  plane  about  the  z-asis.  so 
that  the  x-axis  indicates  north.  This  problem  is  discussed  in  Part  1.  po.  43  ff.  and 
83  ff.  The  basic  stabilization  method  will  be  reviewed  briefly  hero,  after  which 
the  problem  wiil  lie  considered  in  detail,  since  the  self-erecting  latitude  indicator 
to  be  discussed  below,  operates  in  conjunction  with  the  azimuth  stabilizer. 
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The  V-gyro  unit  in  the  vertical  indicator,  that  is,  the  one  w-hos«  input  axis 
coincides  with  the  y  axis,  senses  two  .mpul.tr  velocity  component:  :  the  y-axis 
component  of  the  angular  velocity  ol  tht  indicated  vertical  (normal  to  the  indicated 
meridian  plane),  and  the  component  of  Earth-rate  projected  on  the  y-axis  (see 
hip.  11).  The  y-axis  Earth- late  component  is  a  measure  of  azimuth  misalignment, 
i.e.,  the  departure  ol  the  indicated  east  axis  (the  y-axis)  lrom  true  east  (the  Y-axisj. 
Specifically,  this  angular  velocity  component  is  proportional  to  the  sine  of  the  cer  - 
recti  on  to  indicated  north,  defined  as  the  angle  between  true  north  (the  X-axis)  and 
indicated  north  (the  x-axis).  (Since  the  correction  is  kepi  small  by  the  system,  the 
sine  of  the  correction  is  taken  as  equal  to  the  correction  itsell  in  the  azimuth  stabi¬ 
lizer  synthesis  given  in  Derivation  3.)  The  y-axis  Earth-rate  angular  velocity 
component  can  be  nulled  by  rotating  the  controlled  member  in  the  indicated  hori¬ 
zontal  plane  about  the  z  axis  until  the  y  axis  points  east  and  the  x-axis:  points 
north.  This  is  accomplished  in  the  manner  described  in  Part  1.  Consider  the 
special  case  where  the  angular  velocity  of  the  indicated  vertical  is  zero,  and  an 


Fig.  11.  Horizontal  Earth -rote  component  projection  on  y-ixis  due  to  azimuth  misalignment. 
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azimuth  error  is  present,  with  the  y-uxis  pointed  somewhat  north  (it  .ist.  As 
s  hi  ted,  the  y-gyro  will  sense  the  component  ol  Earth -rate  projected  on  the  y-axis. 

The  y-gyro  unit  then  sends  a  signal  to  the  controlled  meinlier  drive  which  rotates 
the  controlled  member  about  the  >  -axis,  i.e.,  causes  the  x-axis  to  dip  This  in 
turn  will  cause  the  y-axis  accelerometer  unit  to  send  an  orientational  control 
signal  to  both  the  y-  and  z  gyro  units.  The  signal  to  the  y-gyro  and  its  associated 
drive  will  tend  to  level  the1  controlled  member,  while  the  signal  to  the  /.-gyro  unit 
will  pmluce  a  simultaneous  rotation  of  tin  controlled  member  about  the  indicated 
vertical  until  the  angular  velocity  sensed  by  the  y-g.yro  unit  is  substantially  zero, 
with  the  y -axis  pointed  east  (i.e.,  perpendicular  to  the  vector  representing  the 
horizontal  component  of  Earth  rate).  If,  removing  the  restriction  above,  the  base 
has  a  velocity  component  in  a  northward  direction,  the  y-gyro  w  ill  be  nulled  when 
it  is  pointed  slightly  north  of  east  (and  south  of  east  for  southward  velocity).  That 
is,  the  y-axis  will  tend  to  align  itself  in  ine  direciion  of  a  vector  perpendicular  to 
the  vector  resultant  of  the  horizontal  component  of  the  angular  velocity  of  the 
Earth  (in  the  direction  of  Lrue  north)  and  the  northward  angular  velocity  of  the 
base  (in  the  direction  of  true  west).  The  steady -state  indication  of  east  (and  north; 
will  be  in  error  by  an  angle  whose  sine  is  the  ratio  of  the  northward  angular  velo¬ 
city  of  the  base  to  the  horizontal  component  of  Earth  rate.  This  error  will  be  dis¬ 
cussed  later  in  more  detail. 

An  additional  angular  velocity  component  that  is  sensed  by  the  z-gyro  unit  is 
the  vertical  component  of  Earth  rate.  To  prevent  the  controlled  member  from 
being  rotated  about  the  z-axis  at  this  angular  velocity,  an  Earth-rate  compensation 
signal  must  be  provided.  The  effect  of  the  compensation  is  to  maintain  the  gyro- 
unit  axes  in  a  coordinate  system  fixed  w’ith  respect  to  the  Earth  rather  than  with 
respect  to  inertial  space.  The  compensation  signal  which  processes  the  z-gyro  is 
proportional  to  Earth-rate  multiplied  by  the  sine  of  true  latitude.  In  the  practical 
case,  as  given  in  Derivation  3,  this  compensation  signal  is  derived  from  the  \- 
orientational  control  signal  which  contains  Earth  rate  multiplied  by  the  cosine  of 
true  latitude.  This  signal  is  multiplied  by  the  tangent  of  indicated  latitude,  derived 
from  the  latitude  indicator,  to  furnish  the  z-gyro  unit  Earth-rate  compensation. 

The  y-gyro  unit  is  thus  the  signal  source  for  t he*  azimuth  drive  as  well  as  one 
of  two  signal  sources  for  the  y-vcrtieal  indicator  drive,  the  second  source  ir,  the 
latter  case  being  the  y-aeceleronieter.  In  this  sense  the  y-gyro  unit  may  be  refer¬ 
red  to  as  the  east-seeking  gyro  unit,  since,  with  t ho  aid  of  the  v -accelerometer  and 
y-  and  /-controlled  member  drives,  its  irput  axis  tracks  true  east. 

As  indicated  above,  the  y-  and  /-systems  are  coupled,  and  their  dynamic  charac¬ 
teristics  are  interrelated.  Part  I  which  treats  undamped  systems  exclusively, 
shows  (pp.  83  ft.)  the  nature  of  the  interdependence  when  the  /.-system  is  a  second- 
order  system  (i.e.,  representable  by  a  second -order  differential  performance 
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equation.  When  the  y -vertical  indicator  is  damped,  in  the  manner  shown  in  Fig.  2, 
the  coupling  becomes  mure  complex.  It  has  been  found  lea.sible,  in  the  latter  case, 
to  represent  the  /-system  by  a  first-order  equation,  Kq .  (3-211,  in  Derivation  3. 

Tilts  equation  is  applicable  to  the  /.-system  shown  in  Fig.  12. 

The  characteristic  time  of  the  /-system  is  an  important  adjustable  parameter. 
Stability  considerations*  place  a  lower  limit  on  its  value  when  the  y-  and  z- systems 
are  combined  into  a  single  fourth-order  system.**  Figure  13  shows  a  series  of 
transient-response  curves  for  the  coupled  y-  and  /-systems,  Liken  with  the  REAC. 
The  following  assumptions  were  used  in  simulating  the  problem  on  the  analogue 
computer: 

1)  The  y-systeni  is  Schuler -tuned  and  has  a  damping  ratio  of  0.3,  when  de¬ 
coupled  from  the  /-system.  (It  is  assumed,  incidentally,  that  in  the  practi¬ 
cal  case  the  x -system  would  have  similar  characteristics.) 

2)  Component  uncertainties  and  input  inaccuracies  are  ignored,  since  they  do 
not  primarily  affect  stability. 

3)  The  system  is  supposed  to  be  on  a  stationary  base,  so  that  the  oscillations 
shown  in  Fig.  13  occur  about  the  fixed  directions  of  the  true  vertical  and 
true  north,  respectively. 

It  will  be  seen  from  Fig.  13  that  the  general  effect  of  decreasing  the  character¬ 
istic  time  is  to  decrease  the  system  damping.  The  general  effect  of  increasing 
the  characteristic  time  is  to  increase  the  solution  time.  The  entire  system  re¬ 
sponse  is  in  addition  a  function  of  the  true  latitude.  Neither  latitude  variations 
from  zero  to  80  degrees,  nor  /.-channel  characteristic -time  variations  of  compar¬ 
able  effect, have  much  influence  on  over-all  stability  for  a  characteristic  time  of 
about  1.3  hours;  and  this  value  allows  a  solution  time  in  azimuth  of  about  four  hours, 
or  three  Schuler  periods.  The  value  1.3  hours  is  therefore  chosen  tentatively  as 
an  optimum  /-channel  characteristic  time. 

The  effects  of  errors  due  to  imiicrfect  components  and  external  inputs  to  the 
system  deserve  examination.  Equation  (3-25).  the  last  equation  in  Derivation  3, 
displays  these  sources  of  error  when  the  y-system,  /-system,  and  the  self-erect- 
ing-latitude-system  (to  be  discussed  in  the  pages  following)  are  regarded  as  a 
single  system  having  a  fifth-order  performance  equation  in  time -derivatives.  As 
far  as  errors  are  concerned,  (hose  associated  with  the  steady-state  behavior  of 
(he  system  arc  of  greatest  significance,  since  the  other  errors  arc  zero  except 
under  transient  conditions.  There  are  two  such  sources  of  steady-state  error; 


*  Specifically,  if  Roulh's  stability  criterion  is  applied  lo  ‘lie  y-  and  /•combi nrd  systems.  and  the  •. -system  is  ScliuleMuneo  and  bss 
a  damping  ratio  of  0.3  when  decoupled  from  the  .’-system,  the  i.haiarteriMic  tune  n!  'he  /-system  has  a  lewc-t  limit  el  apptoyimatr'v  0.4? 
bout. 

**  the  x-system  will  no!  onto:  Ihese  considerations,  because  of  the  aforementioned  tone, lull  pendulum  study. 
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Fig.  13.  Tronsient  response  of  z-system  (ozimuth)  ond  y-verticol  indicotor  for  successive  z  ond  y  step-inputs  with 
o  ronqe  of  ozimuth  system  chorocteri Stic  time  settings  from  1.9  to  0.54  hour  ot  lotitudes  0,45,  ond  80  degrees. 


1)  A  term  proportional  to  tin1  tn.HHur.uy  tn  the  latitude- rati-  (wait  i  -speed! 
compensation.  This  term  varies  as  tin-  secant  of  latitude.  At  latitude 
•55  decrees  an  iiuceujaey  in  the  latitude-rate  compensation  of  u,7  knot  will 
create  an  error  of  one  mil*  in  indicated  azimuth. 

21  A  term  proportional  to  y-gyro  unit  driit.  Thin  term  also  vane:  a;-  the  m- 
cam  ot  ia 1 1 Hide .  A!  latitude  45  degree;  ,  a  drift  of  0.01  depice  pit  hour  will 
create  an  error  of  one  mil  in  indicated  azimuth. 

The  first  of  these  two  errors  would  involve  the  aforementioned  ratio  of  tin-  latitude 
rate  t  >  the  horrs.iitai  component  of  Karth  rate.  The  u  i  of  a  northerly  water-. ‘■peed 
compensation  signal  reduces  this  error  to  the  residual  one  riven  above-,  which  is 
due  to  the  inaccuracy  in  the  water-speed  compensation.  Note  that  no  cumulauvt 
errors  are  incurred  in  azimuth  indication. 

°0$i  T ION  INDiCA TkON  Br  G-E. GMr  TRIG  iNT t  GRATlON  US:NG  A  Pftf  At.  II. !  [>  inf  r  T  A . 

GYRO  UNiT 

The  position  indication  system  is  shown  in  1-ip.  14,  and  is  accursed  in  Part  I. 
pp.  38  If.  and  101  if.  Tin*  major  objection  to  this  system  i  the  requirement  for 
proalignment  of  the  pyro  axes  with  tli  directions  of  astronomical  :  tar -lines,  dau 
which  are  not  readily  available  with  sufficient  accuracy  on  a  system  in;  hilled  in  a 
submarine.  In  particular,  the  required  azimuth  reference  would  be  difficult  to  (  sta- 
blish,  a  situation  to  be  contrasted  with,  for  example,  the  problem  of  aircraft  guid¬ 
ance,  where  the  reference  can  be  established  on  the  pround  before  take-off.  In 
submarine  navigation  systems,  in  general,  self-erecting  systems  aic  to  be  re¬ 
garded  as  preferable  to  pre -aligned  systems.  Therefore,  prealigned  systems 
are  not  discussed  further  here.  However,  it  will  be  shown  subsequently  that  if  a 
self -erecting  system  is  us  *d  to  orient  an  inertial  gyro  unit  so  that  it  is  s<  nsitive 
only  to  angular  velocity  components  effectively  parallel  to  the  Karth  polar  axis, 
the  gyro  unit  will  be  a  source  of  data  on  longitude. 

A  SFt.  E-ERECTING  GEOMETRIC  LATITUDE  INDICATOR 

The  indication  of  latitude  by  geometric,  integration  was  discussed  in  Part  1. 
pp.  37  ff.  The  integration  is  accomplished  by  mounting  a  gyro  unit  on  a  latitude 
gimhal, which  causes  a  line  on  the  gimbal  to  track  a  line  parallel  to  the  Earth  polar 
axis,  thereby  furnishing  a  reference  vertu  al;  i.i  .,  the  perpendicular  tn  the  Karth 
polar  axis  is  a  vertical  at  the  equator. 

The  indicated  latitude  is  required  in  the  system  itself,  e.g.,  for  vector  com¬ 
ponent  resolution  of  Karth  rate  compensation  in  the  x  system,  for  resolution  of 
the  x- current  when  used  for  compensation  in  the  z-ch  mu  1.  and  for  sensitivity 

*  Tht  mi/  will  fa  fafmfa  *n-  iisir.  report  i‘,  r.  nv»*  f. :  firo-ron|rrj  app1)'#1?!)'*  ,*•  iwp  t  • 
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control  through  the  secant  of  lauiud.  01  the  z-channti.  Thi:  :  ug.gr:-  ts  an  inevitable 
feedback  from  the  latitude  gimbul  to  the  vertical  and  azimuth  lndu  atom.  Tins  feed¬ 
back  ettect  mav  be  utilized  tn  erect  the  latitude  glmbai  us  is  sub:  i-quemly  shown, 
thereby  making  the  latitude  gyro  unit  mentioned  m  Part  !  uniiece  -ary.  The  mech¬ 
anization  of  such  a  device  will  n.;v.  be  described.  The  latitude  system  performance 
equation,  Eq.  (4-1 1>,  is  derived  m  Derivation  4.  Tin  essential  operating  compon¬ 
ents  of  the  latitude  indicator  are  shown  in  Tig.  la. 

The  kind  of  signal  input  to  the  latitude  gimhal  drive  is  determinable  from  the 
principal  requirement  for  indicated  latitude  imposed  b\  tin  re  t  of  the  system: 
namely,  the  multiplication  of  the  x-orientationa! -eontrol  signal  |c>  tin  tangent  of 
indicated  latitude,  so  that  the  product  can  be  used  to  compensate  the  /.-gyro  unit 
for  the  vertical  component  of  Earth-rate  plus  longitude  ehaiipt  .  As  indicated  m 
Derivation  4,  the  latitude  drive  signal  eho  en  consists  of  the  y-  orientational -con¬ 
trol  signal,  with  true  latitude  rate  partly  compensated  by  water  speed,  and.  in 
addition,  the  water-speed  compensation  itself. 

Consider  firs.!  the  case  where  the  base  is  stationary  with  respect  to  the  Kurils 
and  assume  that,  originally  no  errors  are  present  in  the  v<  meal  and  azimuth  sys¬ 
tems,  Then  the  input  axis  of  the  z-gvro  unit  is  vertical  and  the  input  axis  o!  the 
y-gyro  unit  points  east.  Assume  that  the  latitude  gimhal  is  not  positioned  cor¬ 
rectly.  Under  these  conditions  the  latitude  resolver  (.  ec  fig.  16)  will  not  receive 
the  correct  angle  as  an  input.  This  in  turn  will  produce  an  incorrect  resolution  of 
the  x-gyro  unit  current,  which  is  used  to  compensate  the  z-g.vru  unit  for  the  verti¬ 
cal  component  of  Earth  rate.  The  /-gyro  unit  will  then  sense  tin  uncompensated 
portion  of  vertical  Earth  rate,  and,  through  the  azimuth  drive,  will  produce  a  rota¬ 
tion  of  the  controlled  member  about  the  indicated  vertical.  When  this  rotation  oc¬ 
curs,  the  input  axis  of  the  y-gyro  unit  will  be  moved  from  its  eastward  position, 
and  the  unit  will  sense  a  component  of  horizontal  Earth  rate,  f  rom  this  point  the 
action  is  similar  to  that  described  previously  in  t  onneciioii  with  azimuth  stabili¬ 
zation  except  that  there  are  now  involved  three  coupled  systems  instead  of  two. 

The  y-gyro  unit,  by  means  of  its  drive,  will  produce  a  rotation  o|  the  controlled 
member  about  the  y-axis.  This  in  turn  will  cause  the  y  axi:  pendulum  to  send 
simultaneous  orientational  control  signals  to  the  y-gyro.  the  z-gvro  and  the  latitude 
drive.  The  signal  to  the  y-gyro  unit  will  tend  to  level  the  controlled  member,  the 
signal  to  the  /  -gyro  unit  will  produce  a  rotation  of  the  controlled  member  about  the 
indicated  vertical  until  the  angular  velocity  sensed  by  the  y-gvro  unit  is  substan¬ 
tially  zero,  and  the  signal  to  the  latitude  drive  will  i  otuto  tin  latitude  gimhal  (sec 
fig.  16)  until  the  uncompensated  angular  velocity  sensed  by  the  z-gvru  unit  is  zero. 
The  last  condition  will  exist  only  when  the  latitude  gimhal  is  correctly  positioned 
and  tin-  latitude  resolver  receives  the  particular  input  angle  necessary  to  produce 
accurate  compensation  for  vertical  Earth  rate  in  the  /  -gyro  unit. 
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The  operation  of  tin*  ev  ten»  >»n  ;i  I  mm*  with  a  non-zt-i  o  latitude  rate  i.L  or  sen 
tiallv  the  same  except  tor  i!,e  otie  t:  on  the  dynamic  p«  i  im  matte  e  ot  the  azimuth 
stabil i/.atiuti  system  (noted  previously)  ami  the  latitude  :  y.  tern.  'I he  dynamic:  of 
the  latTude  system  untie t  discussion,  when  it  i;  not  rumpem.utod  lot  latitude  l  ate, 
are  such  as  to  make  the  system  •  ubject  to  error  tl  there  t:  a  liniti-  velocity  input . 

1  he  disadvantages  of  such  an  arrangement  are  apparent.  1  w  o  po:  sible  :  oiutrms 
appe.tr  tube  leastble:  vi/.,  to  change  the  system  dynamic:  and  then  by  ptoouce  a 
latitude  system  with  /.ero  velocity  error,  or  to  use  an  external  means  lor  latitude 
rate  eonipeus.i.uw  so  that  the  latitude  g  initial  ts  driven  at  the  approximate  l.ilitudt 
rate  at  all  times.  The  first  methixi  is  superior  tor  the  submarine  pi  iblem.  situ  e 
it  minimizes  the  effect  cf  eorapensation  inaccuracy  and  1  tin  one  which  oiler 
most  promise  for  further  development.  However,  it  does  introduce  a  lability 
problem  which,  because  ol  the  coupling  between  tht  latitude,  azomuth,  and  y -system 
is  difficult  to  solve.  Tilts  problem  is  being  studied.  lor  the  put  pom  >t  the  present 
report,  it  will  be  shown  to  what  extent  a  system  utilizing  water peed  compensation 
is  practical.  Derivation  4  is  therelore  based  on  the  latter  approach. 

In  accordance  with  the  foregoing  discussion,  the  lalitub  drive  current  ts  in  td< 
up  of  the  following  component  signals  (see  Kq.  (4- It): 

])  The  azimuth  orientational  conn  ol  signal,  minus  Earth-  rate  compe-m  atmn. 
This  signal  consists  of  the  y -orientational  control  s*gt...l  plus  latitude- rate 
compensation. 

21  Tn’  rate  of  change  of  latitude,  based  on  pilomotor-log  information. 

Thor  qnals  are  sent  to  the  latitude  gimbal  drive  system,  mounted  on  tht 
controlled  member.  Equation  (4-11),  tne  performance  oquatt  m  lor  the  latitude 
loop,  is  of  first  order  in  time  derivatives  of  the  correction  to  indicated  latitude  ; 
the  equation  for  the  three  systems  menu. •rod,  taken  as  a  single  system,  is  there¬ 
fore  a  quintic  (Eq.  (4-12)).  Routh's  criterion,  when  applied  b  lite  latter  equation, 
with  y-  and  z -parameters  chosen  as  already  indicated,  shows  that  lor  stability  tht 
lower  limit  for  the  characteristic  time  in  the  latitude  system  is  about  0.6  hour. 
Figure  17  shows  tne  response  to  stop -inputs  in  >•-.  z-.  and  latitude-channels  lor  a 
simulated*  system,  it  is  assumed  in  the  simulation  that  the  system  base  is  stall  >n- 
ary  with  respect  to  tin*  Earth  and  that  the  inaccurac  ies  and  urn  ertainties  in  the 
quantities  upon  which  the  system  operates  are  zero;  this  simplification  ts  con¬ 
sistent  with  stability  determination  as  an  isolated  problem.  A  practical  lower  limit 
on  the  latitude-characteristic  time  is  seen  from  Fig.  17  to  be  about  1.3  hours,  and 
this  is  the  recommended  value. 

The  latitude  system  is  subject  to  U  ,  ly -Mate  <  rmrs;  this  is  shown  bv  Kq. 
■,4-i?)  when  the  time  derivatives  of  correction  terms  ate  equated  to  zero.  The 
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input*  and  latitude-system  characteristic  time  settings. 


magnitude  i if  these  eiror:  tan  he  e:  timuted  !m  typical  operating  conditions  by  con¬ 
sidering  the  .submarine  to  be  moving  .it  a  speed  ol  20  kn>>ts  on  a  meridian  course 
.it  latitude  40  decrees.  The  tt  idy-.st.iie  rrmr  in  indicated  latitude  under  the  se¬ 
cond  i  lions  i  made  up  principally  ut  the  following  terms: 

1)  Resolver  uncertainty:  about  1.5  minute,--  o)  arc  error  in  indicated  latitude, 
if  the  resolve]  i.  accurate  lo  i.UO  percent. 

2)  (:yro  unit  drill:  about  1.0  minutes  of  are  error  in  indicated  latitude  for 
either  x-gyro  unit  drift  or  z-gyro  unit  drift  of  0.01  degree  per  hour, 

2)  l.atitude  rate  compensation  inaccuracy:  about  1.3  minute:-  ul  arc  error  in 
indicated  latitude  for  each  knot  ol  inaccuracy  in  latitude  rate  compensation. 

4>  Accelerometer  uncertainty:  about  0.25  minute  of  arc  erm r  m  indicated 
latitude. 

5)  Orientational  control  sipnal  uncertainty:  about  0.2  minute  ol  ai  ■  error  m 
indicated  latitude  for  a  noise  level  in  the  /, -orientational  confioi  ,  ipnai 
generator  ol  0.3  percent. 

Generally  speaking,  in  terms-  of  the  perlorniance  ol  these  components  as  they 
are  used  in  current  applications,  the  principal  serious  sources  of  error  c..n  be-  ex¬ 
pected  to  be  the  resolver,  the  latitude  rate  compensation,  and  the  gyro  units,  ns 
that  order. 

LONGITUDE  INDICATION  BY  CELT  ST  I AL -LONGITUDE  .RATE  TRACKING 

The  indication  of  longitude  necessarily  involves  the  mtepration  ot  longitude 
rate.  A  straightforward  method  of  accomplishiiip.  this  involves  the  use  of  a  pyro 
unit,  coupled  to  a  servo  drive,  as  the  integrator;  tins  at  once  minimize*-  the  effect 
of  integrator  drift  (in  this  case,  drift  effects  are  confined  to  the  gyn  and  allows 
the  use  of  the  latitude  gimhal  just  described  to  orient  the  longitude  gyro  input  axis. 
The  mounting  of  the  longitude  gyro  on  the  latitude  gimhal  is  shown  schematically  m 
Fig.  15.  The  input  axis  of  the  gyro  unit  is  oriented  by  the  latitude  gimhal  in  such 
a  way  that  the  gyro  unit  is  sensitive  only  to  angular  velocity  component.*'  that  are 
effectively  parallel  to  the  Karth  polar  axis  (see  Part  I,  pp.  3H  and  fi?  ff.i, 

These  components  include  the  Karth  rate,  and  since  the  gyro  unit  is  sensitive 
to  motion  with  respect  to  inertial  space,  the  signal  from  the  gyro  unit  to  the  longi¬ 
tude  gimhal  drive  will  cause  the  drive  to  operate  .t  the  c<  lest  in  1  longitude  rate.  A 
separate  time  drive  can  be  used  to  subtract  Karth  rate  from  this  motion,  so  that 
the  longitude  gimhal  rotates  at  the  indicated  longitude  rate,  and  the  angle  between 
the  longitude  gimhal  and  a  preset  reference  direction  on  the  latitude  gimhal  is  in¬ 
dicated  longitude.  It  appears  preferable  to  mount  the  time  drive  m  a  console  olf 
the  gimhal  and  carry  out  the  above  -  ubtractinn  in  the  console.  This  permits  a  re- 
ducti  ,n  in  the  size  of  the  gimhal  and  more  convenient  system  pac  kaging. 
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Tilt-  f  iihitt  mg  conclusions 


Tilth  metlujd  tv  discussed  lit  detail  th  Derivation  b, 
may  be  drawn  front  Derivati on  5,  Kq.  lb-24 

1;  The  transient  behavior  ol  the  longitude  :-.y;  tem  e  .  tmewhai  influenced  In 
the  rest  ul  the  indication  system,  particularly  the  -  and  /  -  systems  .  Him 
ever,  there  is  no  teed back  Iruin  the  longitude  loop  to  the  let  t  ■  <!  tiie  -\st«-n 

2  The  steady -Man  accuracy  nt  tin  system  e  dependent  in  Die  h  mg  tlude  veh 
city,  unless  tiie  longitude  system  rluracten:  tn  nine  is  .*  ullu  n  ml >  small. 

3  i  Tiie  longitude  gyro  unit  contributes  a  drtlt  term  to  the  performance  rqy.i 

non  (Kq.  (5-241.  The  integral  ol  tin: -  unit  term  with  n  pert  to.  time  ap¬ 
pears;  this  means  that  a  eumuiatne  error  can  be  expec  ted,  deb  mimed  bv 
the  drift  rate  and  tune  ot  opet. atoii  nl  tin  sv:  tens. 

•i;  To  a  lesser  extent ,  the  inaccuracy  in  the  -odereal  turn  drne  contributes  !• 
the  above  cumulative  (  rr«»r  tn  the  same  manner  a  gyro  unit  drift. 

5)  The  error  in  setting  reference  longitude  on  the  indicator  dial  is  not  cumu¬ 
lative  anti  appears  in  the  same  form  in  each  longitude  leading, 

.  assess  the  eliect  of  these  error:  on  longitude  indication,  eons  trier  Kq.  ifs-2 
in  Derivation  5  under  steady-state  condition.-  ,  t.e.,  alter  transients  have  settled 
out,  so  that  the  time  derivatives  of  cur  rectum  t»  tons  are  /.ero.  As:  unie  that  th< 
characteri.  tic  time  tor  the  longitude  system  ts  1  n  .•  ecoiKis,  and  that  the  submarine 
is  at  latitude  45  degrees,  traveling  northward  at  30  knots.  Under  the:  e  conditions 
a  one-mil  error  in  azimuth,  or  a  one-mil  error  in  the  x-mmponcnt  of  the  correc- 
turn  to  the  indicated  vertical,  will  produce  approximately  0.02  minnti  of  arc  error 
tn  indicated  longitude  for  each  hour  ot  operatitii’.  time.  A  longitude  gyro  unit  drift 
of  0.015  decree  per  hour  will  be  a  source  ■  f  about  one  minute  of  arc  error  tn  tn- 
viicated  longitude  for  each  hour  of  ope  rat  ini;  time.  If  the  sidereal  time  drive  t«  m 
err-'r  by  0.01  percent  there  will  be  an  error  Iront  this  source  of  about  0.1  minute 
of  arc  in  indicated  longitude  for  each  hour  of  oporatiii|;  time, 

As  stated  initially,  the  indication  of  longitude  requires  that  longitude  rate  be 
integrated.  The  consequence  appears  to  be  an  unavoidable  cumulative  error  as¬ 
sociated  with  the  integration.  The  inevitability  of  the  integration  requirement 
stems  ultimately  from  the  non-unique  quality  of  longitude  and  the  :  ymmetry  of  the 
Earth  about  its  polar  axis,  Tia-refoit  ,  the  engineering  procedure  in  longitude  in¬ 
dicator  design  must  be  directed  at  the  reduc  tion  of  cumulative  error,  .  through  ef¬ 
ficient  disposition  of  components  and  sufficiently  precise  instrumentation, 
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GKOUNO  iPEED  INDICATION 


All  uf  the  positnm-indicatt  ig  .ystems  it*  cnk-d  have  output:  in  tin  lor m  of 
angula r ,  rather  than  iin*ai  ,  tii  .jjS.it  vmcnts  of  tie  submarine  It  is  to  he  *  xpes  ted, 
therefore,  that  ground  i  obtained  from  tin.  type  ut  navigation  system  will  be 
derived  from  the  essentially  geocentric  angular  velocity  of  the  lndieatt  u  vertical 
with  respect  to  the  harth  The  product  of  the  quantity  aim  the  lauiu.  of  the 
Earth  then  gives  the  ground  ;  peed.  '1  hut  the  chief  problem  in  indicating  ground 
speed  is  the  extraction  of  the  angular  velocity  of  tin  indicated  vertical  with  it 
spect  to  the  Earth  front  tin  navigation  y  t<  m 

This  angular  velocity  may  be  obtained  either  trotn  11  the  vertical  indicator 
itself  or  2)  the  latitude  and  h  -igituue  indicator:  fhe.se  methods  an  dcsci  lx  o 
in  Derivation  6.  Ground-speed  indication  by  each  of  these  two  methods  will  now 
be  disctiss..<l 

IMF.  VERTICAL  INDICATOR  as  a  source  of  ground-spec d  data 

hi  method  ill  where  the  vertical  indicator  is  the  initial  source  of  ground - 
speed  data,  the  quantities  ol  importance  are  tin  x-  and  y  axis  orientational  con¬ 
trol  signals.  These  signals  an-  respectively  proportional  to  the  x-  and  y- com¬ 
ponents  of  the  angular  velocity  ol  the  indicated  vertical  with  respect  to  inertial 
space  The  y  -  component  of  the  angular  velocity  of  the  indicated  vertical  wit!;  r<  - 
sped  to  the  Earth  is  equal  to  that  with  respect  to  inertial  space  except  for  small 
correction  terms.  It  is  therefore  at  once  u  aide  to  indicate  the  north*  *  1\  ground- 
speed  component.  The  x-eomponent,  however,  contains  an  Earth- rat*'  term, 
which  must  be  removed  be  for «  east-west  ground  speed  is  available  from  this 
source.  Specifically,  a  signal  proportional  to  Earth- rate  compensation  multiplied 
by  tiie  cosine  of  indicated  latitude  (the  latter  derived  from  the  latitude  indicator! 
is  subtracted  from  the  x-axr-  orientational  control  signal,  i  lie  resultant  signal 
i  -  approximately  proportional  to  the  x  component  of  the  angular  velocity  of  the 
indicated  vertical  with  r*  meet  to  tin  Earth  The  indicated-latltudi  function  and 
Earth-rate  compensation  loin  cause  the  \ -component  to  be-  influenced  to  som* 
extent  by  Earth  rat*-  compensation  inaccuracy  and  latitude-  rate  compensation  in¬ 
accuracy.  Similarly,  the  v  component  **f  the  angular  velocity  of  the  indicated 
vertical  with  respect  to  th*  Earth  is  affected  by  the  correction  to  indicated  north 
These  errors  arc  ubsequentlv  discussed;  for  tin  present,  note  that  two  signal 
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can  lie  divided  from  the  ventral  indicator,  each  ess,i  ntiailv  proportional  to  a 
component  of  ground  speed. 

It  is  necessary  that  these  components,  be  added  m  quadrature  tiv  a  resolver 
set  at  an  angle  o!  90  degrees,  so  that  the  r«  solve!  output  will  tit'  pi  opui  tional  to 
the  magnitude  of  the  vector  angular  velocity  of  the  indicated  vertical  vtilti  respeci 
to  the  F  arth.  The  result  of  such  addition  is  given  for  Uie  ca.  e  of  a  system  lu 
equilibrium  by  F.‘q  (6-131  of  Derivation  ft.  The  indit  ated  ground  speed  ss  equal  to 
the  true  ground  speed  plus  an  inaccuracy  which  depends  on  the  direction  of  mo¬ 
tion  over  the  F  arth,  1  e.,  on  the  cosine  or  suit  of  true  azimuth.  Tin  impoiiant 
influences  on  the  inaccuracy  term  are  due  to: 

1)  The  x-gvro  and  y-gyro  drift  terms:  lor  each  0.01  deg  ret  per  horn  of 
drift,  the  contributory  error  in  the  eastward  id  northward  indicated 
ground-speed  components  respectively  will  he  approximately  0.6  knot 

21  The  term  containing  the  correction  to  indicated  north;  at  latitude  45  de¬ 
grees,  each  rml  ot  error  in  indicated  azimuth  contributes  approximately 
0.6  knot  of  error  to  the  eastward  indicated  ground-speed  component  Tin 
contribution  to  the  northward  ground-speed  component  error  from  tins 
source  is  much  smaller,  being  only  about  0.02  knot  for  each  mil  of  izi- 
mutli  error,  at  a  ship  speed  of  20  knots.) 

3l  1  he  Harth-  rate  compensation  inaccuracy  term;  it  the  compensation  •  lgnal 
is  constant  to  within  an  accuracy  of  0  01  pern  nt,  the  contributory  <  m>r 
m  the  eastward  indicated  ground-speed  component  will  lie  about  0.09  knot 

41  The  term  containing  the  correction  to  indicated  latitude:  latitude*  45  d<  - 

grees,  each  minute  of  arc  error  in  latitude  indication  contributes  about 
0.2  knot  of  error  to  the  northward  component  of  indicate  d  ground  :  pci  d 
The  y-axis  correction  to  the  indicated  vertical  contributes  the  same  rela¬ 
tive  error  to  this  ground-speed  component  and  the  x-axis  correction  to 
the  indicated  vertical  contributes  the  same  relative  t  m.u  to  the  ea:  twarri 
component  of  indicated  ground  pi-ed. 

THE  LATITUDE  AND  LONGITUDE  INDICATORS  AS  SOURCES  Of  GROUND*  SPf  l  D  DATA 

If  ground-speed  information  is  obtained  from  the  latitude  and  longitude  nidi 
r  a  tors  f  method  i2)  above),  it  is  necessary  to  add.  m  quadratun  .  signals  propor¬ 
tional  to  the  indicated  longitude  rate,  multiplied  by  the  cosine  of  indicated  latitude . 
and  the  indicated  latitude  rate  One  of  the  signal.*  -uitahle  for  thi,  purpose  is 
that  furnished  to  the  latitude-  drive;  tins,  signal  is  proportional  to  indicated  lati 
tude  rate.  The  indicated  longitude  can  be  obtained  from  a  tachometer  attached  to 
a  shaft  whose  rotation  rate  is  proportional  to  indicated  longitude  rati  The  indi¬ 
cated  latitude  is  obtained  as  before  from  the  latitude  system.  The  indicated  ground 
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.speed  is  then  given  by  i.q  »6-20i  ul  Ijci  nation  '  I  In  angular  velocity  »1  the 
Indicated  vertical  is  equal  to  the  angular  velocity  of  the  ‘rue  vertical  plus  an  in¬ 
accuracy  term  as  be!or«  the  luipui  lance  influences  on  the  inaccuracy  in  this 
case  being  as  toiiows: 

1)  I  he  longitude  gyro  drill  term,  at  latitude  45  degrees,  lor  cadi  0  01  de¬ 
gree  per  hour  ul  drift,  the  contributory  erroi  in  the  northward  gtoutid- 
speed  coinpoiient  will  hi  appi  oxunateiy  0  4  knot. 

2,i  The  longitude-rate  tachometer  uncertainty  term:  this  contributes  an  error 
to  the  northward  indicated  ground-speed  component ,  at  latitude  4f>  degrees 
and  with  a  component  of  20  knots  magnitude,  ol  about  0.1  knot  tor  a  one 
percent  uncertainty  in  the  tachometer 

3)  The  l.arth  rate  compensation  'Sidereal  turn  drive)  inaccuracy  term  at 
latitude  45  degrees,  when  the  error  in  the  tiine-drivt  rotational  rate  is 
0.01  percent,  the  contributory  <  trot  in  the  northward  indicated  ground - 
speed  component  will  be  about  o  00  knot  Note  that  the  sidereal  tune-drive 
inaccuracy  appears  in  the  inaccuracy  in  the  indicated  ground  speed,  a 
compared  with  the  Karth-rate  compensation  signal  inaccuracy  m  the  pr«  - 
vtous  case.  At  present,  a  tune-drive  can  be  engineered,  in  general,  as  an 
inherently  more  accurate  device1  than  a  constant-level  signal  generator 

4)  Terms  containing  the  correction  to  indicated  latitude,  the'  corn  ction  to 
indicated  north,  and  the  x-axts  correction  to  the  indicated  vertical;  at 
latitude  45  degrees,  and  northward  and  eastward  ship  velocity  compo¬ 
nents  of  20  knots  each,  one  mil  of  error  iti  azimuth  place;  the  northward 
ground-speed  component  in  error  by  about  0.01  knot,  white  one  minute 

of  arc  error  iri  latitude  or  x- vertical  indication  contribute;  approximately 
0.005  knot  to  the  indicated  ground-speed  error  northward  component  Note 
that  at  equilibrium,  l.c  ,  whe  n  the  systi  m  is  not  suffering  a  transient 
shock,  the  eastward  ground-speed  component  is  without  errors  of  the  type 
discussed  above. 

One  advantage  of  this  method  over  the  first  one  discu:  sod  is.  that  the  gyro 
drift  is  here  multiplied  by  the1  cosine  of  latitude.  Furthermore  .  the  azimuth  e  rror 
which  is  finite  in  the  steady  state,  makes  a  less  serious  contribution  hire  It 
may  be  concluded  that  the  method  <jf  ground-speed  indication  discussed  imme¬ 
diately  above  is  preferable  to  the  first  method.  The  principal  limitation  would 
derive  from  the  requirement  for  a  precise  tachometer  for  the  measurement  ol 
longitude  rate 
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RECOMMENDED  SUBMARINE  INERTIAL  NAVIGATION  SYSTEM 


GENERAL  CONSIDERATIONS 

Heport  H-9,  Parts  1  and  11,  represents  the  result  of  a  study  of  .several 
method  of  vertical-  and  position -unite ation.  of  winch  the  most  promising  have 
been  discussed  The  basic  theory  has  been  presented,  along  with  suggestions  as 
to  hoi*  presently  available  components  could  be  assembled  into  a  workable  .-as  - 
tern 

The  nature  oi  .he  instrumentation  art  as  reflected  in  the  precision  of  compo¬ 
nents  currently  under  development,  rather  than  any  basic  theoretua!  reason,  de¬ 
termines  the  choice  of  the  best  navigation  system  from  the  engineering  point  of 
view  Four  classes  of  components  of  particular  significance  in  their  elfecl  on  the 
performance  of  an  inertial  position  indicating  system  an  : 

1)  gyro  units, 

21  accelerometer  units, 

3)  latitude  resolvers, 

4)  sidereal  time  drives. 

The  present  precision  of  these  units  may  be  specified  as  follows:  a  gyro 
drift  of  0.01  degree  per  hour  is  within  the  goal  of  current  system  development, 
accelerometers,  used  in  a  vertical  indicator  of  the  type  discussed  m  this  report, 
can  be  expected  to  contribute  about  0.25  minute  of  arc  error  in  indication  of  the 
vertical;  a  resolver  accuracy  of  0.05  percent  is  attainable;  and  a  s.oereal  time  drive, 
such  as  is  required  in  a  longitude  indicator  using  a  gyro  unit  as  an  integrator, 
c;  n  be  constructed  as  a  crystal -controlled  clock,  with  a  shalt  rotation  that  is  pro¬ 
portional  to  Earth  rate  to  within  one  part  in  1 00.000. 

It  has  been  shown  in  this  report  that  in  certain  respects,  some  position  in¬ 
dicators  are  superior  to  others.  The  selection  of  the  best  position- indication 
method  necessarily  involves  a  compromise  because  no  single  method  is  entirely 
free  from  faults,  ana  even  the  best  systems  described  here  are  susceptible  to 
improvement  through  further  study 

VERTICAL  INDICATION.  Roll  ond  e,,eh 

As  a  result  of  this  study  it  appears  practical  to  make  each  component  sub¬ 
system  of  the  vertical  indicator  similar  to  that  shown  in  Fig  2.  The  optimum 
system  parameters  are  given  on  page  14.  and  the  theory  is  based  on  Derivation  i 
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tn  the  Appendix.  The  controlled  member,  nr  azimuth  gimbai,  i  t < »  in  mounted 
inside  .i  pitch  gimhal,  and  tin.-.  in  turn  inside  a  roll  ginibu! .  Such  an  a:-:  t  mbiy  i: 
sketched  tn  detail  tn  Mg.  18 

I  he  dynamics  of  this  vertical  indicator  are  adjusted  so  mat  it  i  Si  hulei 
turn  t.  1  e. ,  the  eighty -tout  -  ml IJUle  pellod  ctl.H  ai  ti  l  l:  lie  s  .li  e  !!ii  OI  j.’uJ  ait  d  It: 

|)<- riormance  depends  un  thi  property  lather  than  on  precise  exHinai  compctisu- 
t ion  inputs.  Krrors  -Alii  be  o!  an  csscnttaliv  dynaiiiie  natui  i  ,  amt  do  not  accumu¬ 
late  With  tlllle  Ihe  chief  source  of  i  Mol,  due  to  the  accelerometer,  Is  e.\,.’tited 
Id  produce  about  0.2b  minute  of  an-  error  tn  vertical  indication  alter  ihe  tran¬ 
sient  sh.ge  is  over  1  lie  solution  time  (the  time  required  to  m  Hie  out  t  associated 
111111  the  transient  .stage  is  about  ihi  ei  hours . 

By  means  of  the  gtniba!  system,  the  vertical,  i  «  ,  she  azimuth  gutibal  orienta¬ 
tion,  is  compared  with  the  roll-  ami  p'tch -t  iinlia!  orientatum^  to  imltcati  i<  - 
spectively  the  ruff  alii!  pitch  o!  the  ship 

A  Z  :  MU  7  H  VI  AHiijy  A  T  :0N,  De-.it.lny 

I  he  system  of  Mg.  12  is  recoiiitiiended  for  azimuth  stabilization.  1  tie  funda¬ 
mental  operation  involved  i  the  nul lint*,  of  the  angular  velocity  input  to  the  ea:  t 
seeking  gyro  unit  used  bv  the  Vertical  indicators  on  the  azimuth  gimbal.  Thi 
operation  i  analogotr  to  that  of  a  standard  gyrocompass,  orient i tip  the  azimuth 
qimh.il  so  that  a  reference  line  on  the  i; inil.ial,  perpendicular  to  the  input  axis  of 
the  east -seeking  gyro  unit,  will  Hack  true  north  when  latitude- rate  compensation 
is  provided  from  an  external  ource.  1  lie  inaccuracy  m  thi:  compensation  is  a 
source  of  steady-state  error,  to  the  extent  of  about  oue-mii  error  m  indicated 
azimuth  for  a  0. 7-knot  error  hi  latitude- rate  compensation,  at  latitude  45  degrees. 
The  drift  of  the  east -seeking  gyro  unit  under  similar  condition:  will  it  uit  in  an 
error  of  one  mil  in  indicated  azimuth  for  a  gyro  unit  drift  of  0  01  degree  per  hour. 
‘I'll*’  solution  time  of  the  azimuth  s\  ti  ni  i:  about  four  hour  Tin  optimum  sy: 
teni  parameters  are  given  on  page  28  ,  and  tin  tlieon  tical  consul!  rations  an 
given  in  Derivation  3  Figure  18  show  this  azimuth  system  added  to  the  azimuth 
gimh.it.  The  angle  between  the  north-referenre  Inn  on  the  azimuth  gitnbal  and 
the  fore-aft  line  of  the  ship  is  indicated  ship  heading 

LAHTUDF  INDICATION 

I  he  -elf-erecting  latitude  indic  ator  described  on  page 32  is  recommended. 

The  theory  is  given  m  Derivation  4.  It  involves  the  rotation  of  a  la  solver  whose 
output  is  fed  hack  into  the  azimuth  sv  item  Since  the  azimuth  and  ve  rtical  indi¬ 
cators  are  coupled,  a:  mentioned  above,  latitude  indication  is  ultimately  momt  >red 
bv  the  accelerometer  unit  associated  with  the  east-we  t  axis  on  the  vertical  indt 
r  ater.  Tk  latitude,  azimuth,  and  east-we  t  vertical  ind);  ators  make  up  a 
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Fig.  18-  Pictorial  cutaway  of  gimbal  configuration  for  vertical  and  aiimuth  »y»t*r.i». 
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.  19-  Pictorial  cutaway  o<  latitude  and  longitude  gimbul  configuration. 


i  ioseu-  loop  system.  ana  latitude  indication,  along  w  itii  *.<  final  anti  azimuth  indication, 
1  not  so-  ;•  ,:t  to  cumulative  error;  The  rotation  •>*  tin-  latitude  resolver  ts  as¬ 
sociated  with  the  rotation  of  the  latitude  jdmbai,  wh'cli  ef f*  (  lively  track.-  the 
l.artti  polar  axis,  and  is  sketched  in  Fig  19  Hits  system  ts  such  that  the  latitude 
gimbal  is  seif-erti.'in j,  and  rotate  until  a  reft  r«->  it  line  on  the  gtmbal  indicates 
the  polar  axis  without  initial  setting  or  preatignmei  t  latitude- tatt  eompens  a - 
turn  from  an  externa!  source  ts  required  to  minimise  a  servo  velocity  error 
Indicated  latitude  will  be  in  error  by  approximately  one  minute  of  arc  for  t  ach 
knot  of  error  in  the  i.ttiiudt  -rate  compensation,  and  by  approximately  1.5  mtnutes 
of  arc  for  a  drill  of  0.01  degree  per  hour  ol  tin-  two  gyro  unit:  whose  input  axes 
ate  m  the  indicated  meridian  plane.  A  0.05  percent  uncertainty  m  tin  output  of 
the  latitude  resolver  will  also  create  an  error  of  about  1.5  minutes  of  arc  in  in¬ 
dicated  latitude  These  are  steady-state  errors,  alter  the  transient  stage  has 
passed  The  solution  time  is  atxiui  tour  hours.  The  optimum  latitude  system 
parameters  arc  given  on  page36  The  angle  between  the  indicated  horizontal 
(indicated  by  the  azimuth  gtmbafS  and  the  indicated  polar  ax:.-  (furnished  by  the 
latitude  gtntbal)  gives  the  indicated  latitude. 

LO»Ci!  Tu0£  iNOtCA  T  ION 

The  most  practical  method  for  indicating  longitude  is  to  track  celestial  longi¬ 
tude  rate  with  .hi  integrating  gyro  unit  mounted  to  rotate  within  the  latitude  gtm- 
bal,  as  shown  in  Fig.  20.  This  gyro  unit  re  ponds  to  the  angular  velocity  parallc  1 
to  the  indicated  F  irth  polar  axis.  The  gt tribal  on  which  the  longitude  gyro  unit  is 
mounted,  rotating  about  the  indicated  polar  axis  within  the  latitude  gimbal,  forms 
the  innermost  of  the  five  concentric  gimbals  of  the  navigation  system.  The  gim- 
bai  rotates  at  the  indicated  celestial  longitude  rati  I  in  information  is  earned 
to  a  mechanical  differential  remote  from  the  gimbai  system,  where  Farth  rate, 
generated  as  a  shaft  rotation  by  a  .stderc»i  elm  k,  is  subtracted  from  indicated 
celestial  longitude  rate  Tin  differential  output  gives  indicated  longitude  difft  r- 
■  nee  This  output  is  read  as  indicated  longifude  relative  to  a  longitude  reference 
index  setting  on  the  indicator  The  chief  errors  in  longitude  indication  art  cumu¬ 
lative  Under  typical  operating  conditions,  a  longitude  gyro  unit  drift  of  0.015 
degree  per  hour  will  be  a  source  of  about  one  minute  of  arc  error  in  indicated 
longitude  lor  each  hour  of  operating  time.  The  sidereal  time  drive  may  be  ex¬ 
pected  to  contribute  an  error  in  indicated  longitude  of  about  l.;3  minutes  of  arc  in 
a  week  of  opera  tig  time  Tim  longitude  system  is  discussed  on  page  3R  and  in 
Derivation  5 
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GROUND-SPEED  INDICATION 

It  is  recommended  that  ground  speed  be  computed  by  utilizing  signals  pro¬ 
portional  to  indicated  longitude  rate  (from  a  tachometer  geared  to  the  longitude 
indicator)  and  to  indicated  latitude  rate  (the  latitude  gimbal  drive  input  signal). 
The  indicated  longitude- rate  signal  is  then  multiplied  by  the  cosine  of  the  indi¬ 
cated  latitude,  obtained  from  a  resolver  on  the  latitude  gimbal.  This  quantity  is 
then  added  in  quadrature  to  the  indicated  latitude-rate  signal,  by  means  of  another 
resolver,  to  obtain  the  magnitude  of  the  essentially  geocentric  angular  velocity  of 
the  indicated  vertical.  This  last  quantity,  multiplied  by  the  Earth  radius,  is  the 
indicated  ground  speed.  The  chief  error  in  the  indicated  ground  speed  is  in  the 
northward  component,  to  which  a  0.01  degree  per  hour  longitude  gyro  unit  drift 
will  contribute  an  error  of  about  0.4  knot.  This  method  of  ground-speed  indica¬ 
tion  is  discussed  on  page  41  and  in  Derivation  6. 


GENERAL  DATA 

The  following  tabulation  summarizes  the  over-all  characteristics  of  the  pro¬ 
posed  system: 


Approximate  outer  dimensions  of  gimbal  42  in.  high,  42  in.  wide,  and 
system,  including  cover  46  in.  long 

Radius  of  largest  gimbal  (to  go  through  36  in. 

a  standard  20  in.  by  39  in.  hatch,  during 
installation) 


Number  of  self-contained  single-axis  in¬ 
tegrating  gyro  units  required 

Number  of  accelerometers  required 

Number  of  integrators  required* 

Estimated  input  power  required  from  ship 
nower  supplies 

Estimated  dimensions  of  console  mount¬ 
ing  electronic  components 


4  (3  inertial  gyros,  1  used 
as  an  integrator) 

2 

4 

5  kilowatts 

5  ft  high,  5  ft  long,  and  3  ft 
deep 


•  The  number  of  integrators  might  be  reduced  to  three  by  using  the  /-system  in  its  entirety  as  one  ot  the  y- system  integrators.  This 
line  ot  attack  is  being  investigated. 
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NUMMARY  AMD  RECOMMENDATIONS 


The  submarine  navigation  system  discussed  in  the  preceding  section  and  pic¬ 
tured  in  the  sketches,  Figs.  18  and  19,  is  shown  as  a  wooden  mock-up  in  the 
photograph,  Fig,  20.  While  this  system  is  regarded  at  present  as  most  suitable 
for  the  application,  its  recommendation  is  also  described  as  tentative  because  it 
is  based  on  a  time-limited  study.  Nevertheless,  it  is  possible  to  suggest  the  di¬ 
rection  in  which  subsequent  work  might  best  proceed  in  order  to  improve  on  the 
proposed  system. 

Four  gyro  units  are  required  in  the  proposed  system;  three  of  these  are 
inertial  gyro  units,  used  to  set  up  a  three-axis  coordinate  system,  while  the 
fourth  (the  longitude  gyro  unit)  is  used  as  an  integrator.  The  first  three  gyro 
units  represent  an  irreducible  minimum  number  if  single-axis  gyro  units  are 
used.  It  may  be  possible  to  combine  the  function  of  the  longitude  gyro  unit  with 
that  of  the  other  units,  and  reduce  the  total  number  of  single-axis  gyro  units  to 
three.  This  would  require  a  somewhat  different  gimbal  configuration  from  that  of 
the  system  proposed  in  this  report. 

Another  field  for  future  work  based  on  the  present  study  is  in  the  problem  of 
the  adequate  use  of  external  compensations.  Particular  reference  is  made  to  the 
water-speed  compensation,  a  signal  derivable,  for  example,  from  the  ship  pi- 
tometer-log.  The  eastward  and  northward  components  of  the  water  speed  are 
used  by  the  two  component  systems  of  the  vertical  indicator.  The  specific  appli¬ 
cation  here  is  the  minimization  of  acceleration,  not  velocity,  error;  the  vertical- 
indicating  system,  in  effect,  serves  to  differentiate  external  water-speed  informa¬ 
tion  in  using  it  to  compensate  for  dynamic  errors.  The  result  of  this  particular 
application  of  water-speed  data  is  to  leave  residual  errors  in  the  vertical  indica¬ 
tion  that  are  proportional  to  the  rate  of  change  of  the  inaccuracy  in  the  water- 
speed  compensation  signal  rather  than  to  the  inaccuracy  itself,  i.e.,  no  steady- 
state  velocity  error  stems  from  this  source.  Furthermore,  since  the  compensa¬ 
tion  is  dynamic  in  character,  its  effects  arc  strongly  dependent  upon  frequency, 
only  oscillations  near  the  Schuler  frequency  are  of  consequence,  and  for  frequen¬ 
cies  outside  this  region,  the  compensation  is  not  important. 

In  the  proposed  navigation  system,  these  considerations  do  not  hold  for  the 
azimuth  and  latitude  indicators.  Water-speed  compensation  signals  are  used  in 
these  cases  essentially  to  eliminate  effectively  geocentric  angular  velocity 
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components  of  the  indicated  vertical  in  the  process  of  determining  the  direction 
of  the  Earth-rate  vector.  The  result  of  this  app'ication  of  water-speed  compen¬ 
sation  is  steady-state  errors  in  indicated  azimuth  ana  latitude  that  are  propor¬ 
tional  to  the  inaccuracy  in  the  water-speed  compensation,  not  to  its  rate  of  change. 
It  is  therefore  desirable  to  investigate  other  means  to  remove  this  source  of 
steady-state  error,  and  thereby  relieve  the  present  dependence  of  the  azimuth 
and  latitude  systems  on  data  external  to  the  navigation  system.  A  possible  ap¬ 
proach  is  to  feed  indicated  latitude  rate  from  the  latitude  system  back  into  the 
azimuth  indicator.  This  involves  additional  complexity  in  the  concept  of  the  sys¬ 
tem  as  a  w*  lie,  and  possibly  raises  stability  problems.  This  approach  is  a  sub¬ 
ject  of  prcst.it  studies,  and  appears  promising. 


APPENDIX 


The  mathematical  derivations  on  the  following  pages  are 
referred  to  by  number  in  the  foregoing  text.  The  self-defining 

notation  used  in  the  equations  is  based  on  a  formulation  given 

(71 

by  Draper'  \ 
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DERIVATION  1.  DAMPED  SINGLE-AXIS  VERTICAL  INDICATOR. 

A.  DERIVATION  OF  DIFFERENTIAL  PERFORMANCE  MUTATION. 

If  1(C)  V  ft  i)  *s  as  the  correction  to  the  indicated  vritical, 


Pl^)V  (t,i)  '  *(FV)t  "  *(KV)i 


U-l) 


W(FV)i  ’  ^(cmdsXi.W  > ‘{hfi)  •  fl>Wp)  *  <“>"(*. 


<J-2> 


Refer  to  Fig.  1- 1  • 


Vhfil  *  ^(hCKe/i)  p  l*(dir)  4  e(Ifi)  ’  C(cp)  ~  “(ft')1  4  ^'(hfiXdf)  '  *{hfi)o  l'1_3i 


*(  fir)  =  S<diri(e,e)  e(au)  *  ^U^e(d 


Mir) 


(1-4) 


'(Ifi) 


*  ViXe.e)  *(au)  *  (UKtfiXdf)  4  e(lfi)o 


n-4) 


e(cp)  5  (PF\cpX«,e)  W(FVKcp) 


(1-f-) 


where  W(py)(cp)  ‘s  the  equivalent  angular  velocity  of  the  true  vertical  derived  from  pitomrtcr-log  data. 

“d  (PF>(cPX*,e)  is  to  be  determined. 

effb)  *  S(fbXi.e) '(hfi)  4  *UKfb)  !1"‘) 

e(au)  4  S(auXa,e)  {9  [(C > V!  (t,i)  *  RPW(FV)t  }  4  <U>«(au)  <*-*> 

Substitute  the  foregoing  in  Eq.  (1-2).  The  result  is 

S(hf  x 

^(FV)i  *  ^(cmdsXi,®')  ^  _  t  «  -  «  ^(dirXe.e) ^(auXa.e)  9  (OV  (t  ;) 

P  S(hfi)(e,i)SffbXi,e) 

4  RPW(FV)J  4  S(dirXe.e)(UKau)  *  ^U^e(dir)  4  S(lfiXe,e)^au)(a,e)p  l9  UC) Vj{(_ 


RpW( 


(FV)t 


*  S/if,- 


1 


1 


'(IfiXe.e)  p(U)*(au)*  (UHuiVdfi  4  *rifiv,  4 


% 


(lfiXdf)  e(lfiX>  4  '(cpX«.e)"(EVXcr> 


«••»»>  • . .  s,emrx:r •  «w,  •  <»v 

P  5(hfiXe,i)5(fbXi,e) 


(1-0) 
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flf&f&tHa* 


^cmdsXi.T)^(WiXe,i^^  VKcp) 

^ViKa,«)  1  ....  r  ,  r  ?» i  \ 

"  t  ’P  '  tr  T\  ’'U,e!aui  ‘  ^cmdsXi,»nh«Xe,i)P  l'U*(dw) 

TauKa.e)  VV"  1  fyd) 

~  ^e(fh)  ~AcmdsXi,*)^!MiX*.i)^*UfiXdf> 

5(tmJ*Xi,»‘)P^)*(hfiXdO  ~  ^  *  irj)  ‘^W(1F  K<  p>  * 

1  <>K>  (1-1 

I  hr  jerk  error  is  minimized  through  Schuler  tuning,  i.e.,  the  coefficient  of  P*^j •  \’)t  IS  PGuaterf  lc  ^ero. 

1  here  results 


*(ViXa,i)  15  p 

s 

where  R  is  a  set  value  for  the  larrh  radi 


u-n> 


Assume  ?  R  in  twe  following.  The  acceleration  error  is  minimized,  in  part,  by  using  the  pjtometrr- 

log  data,  as  follows.  Define 

fD'W,)  5  W(FVKcp>  -  WUV)t  "*“> 


(pp>(cr.x*» 


r  l  _  L_., 

.  (GT>(U)  (CT)(id> 

^(cnds)(i,T)^(hfi)(r,i) 


(1-0) 


The  differential  performance  equation  is  then 


'icTil;.pi  ’l’  'icn.,,  *V',c’v- 

l(CT)(lg)  “  (CT)(ldp  p(l)W^v*fP>  "  pI  *(CT)!le)P  (l)WnF.Xcp) 


♦  (U)W,  ,}  - 


(KU’  "R.WJ)  (CT)< 


“  ^cmdsXt,T)^hfiXe,i)P  i^*(dir)  ^e(fb>  : 
'  ^(cmdsXt,H)^(h:iXe,i)^^*(lfiXdf) 


S(cmdsx;,”;o(U>ich.'.xdf) 


u-  l<n 


H.  ( >!' TIMI / ATiON  OF  UNDAMPED  RESONANT  FRI  O!  I  \(  V 


Considering  only  the  water-speed  compensation  rffoi  term  in  f  q,  the  equation  be 

1  Q  Q  1 

{ p J  ♦  p 2  ♦  p*  (C )  V  /  • 

(CT)f1?)  P  R,  P  Rs  <CT)(U) 


‘  (CT)Us)  ’  (CT)(l  i)  p(,)W<l  vKq- 


!>) 


( J  -  D 


The  problem  is  to  select  that  value  for  the  resonant  frequency  associated  with  a  selected  viamphn*  Mtio  of 
the  system  for  which  the  right'hand  side  ot  the  equation  is  a  minimum.  To  do  this,  the  follou sne  notation 
is  defined 

W  12  undamped  resonant  frequency  of  the  quadratic  term 
(DR)  damping  ratio  of  the  quadratic  tern 
(CT)  =  characteristic  time  of  first-order  term 

g 

W  5  =  thr  Schuler  natural  frpqocncv 


1  hen  in  operator  form 

1 

(CT), 


W2  n  ‘ 


W2 

ns 


<1/0 


nsr  (CT)(W)  (CT) 


*  (P  */rTJ  <P*  ♦  ;2(DR)Wn'  P  ♦  VM  ) 


2(DR)Wn 


W2 


=  o’  *  ^2{DR)W  *  — = —  )p*  ♦  *WCp  ♦ 

P  "  (CT)  P  (CT)  n  P  (CT) 


n-isi 


Offine 


(CT)(W  (CT)f|d) 


From  Fq.  (.1-181 


1  1  w 2 

F  «  2(DR)W  * -  -  - 

n  (CT)  (CT)  r 

ns 


n-i<n 


and 


2(0R)Wn 

(CT) 


W2  =  w2 

n  ns 


fl-?01 
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C.  TRANS3’  NT  RESPONSE. 


The  homogeneous  equation  derivable  from  Kq.  (1-1A)  is 

{p3  t  -  1  p2  4  ,9  p  ♦  .9  1  }  f  {C  )  V 1  ,  .  -  Q 

P  (CT)W  P  R.p  R.  (CT)(U/ 


'(id) 


This  is  of  the  form 


2(DR)W  W2 


lpJ  ♦  [2(DR)Wn  *-]p2  *  (W2  t--^]p+_-3 


The  solution  is 


[(C) V](  ;)  *  d(C)V](ti))0{(_ - - 

(MJ  1  2(DR)W 


(CT)2  (CT) 


♦  W2 


1  _ 

+  CT  Vi  “(DR)2  ;  1 

(ctF1 


1 _ 

2(DR)Wn 

(CT) 


x  e  -(»)»,>  sin(WnVl  -  (DR)2  t  *  A(ph)1  j 

2(DR)W„ 

♦  ([(C)V](t>i)I 


°  1(_1_  2(DR)Wn  )e 
(CT)2  (CT) 


1 


Wn  Vl  "  (DR)2 


1  2(DR)W 

- - —  +  +  W2 

(CT)2  (CT) 

_J_.  -  2(  DR)  W  W2 

(CT)2 - -  +  " 

"•  ’  (CT) 


x  e  "(nR)wnt  sin(Wn  _  (DR)2,  ♦  A(ph)2  ] 


(1-27) 


[(C) V)<t  i)  *  0  (1-28) 

"  (CT) 


t 

(CT) 
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where 


MUC)V)(t,0)o  {(-j:  ~T(DR)w7  W2 ) ' 
(CT)2  '  (CT) 

,  1  (  V 

Wn(l  -  (DR) 2  \  1  2(DR)Wn  W2 

(CT)2  "  (CT) 


e  "(DR)V  «nWnVl  -  (DR)2t  -  A(ph)}  } 


♦ 

(CT) 


\phH 


tan 


Vl  “  (DR)2 

TdrT  ton 

wn  VI  -  (DR)"2 


Wn  Vl  (DR)2 


CT  -  (DR)wn 


V)2 


♦an 


«*>»„  'cY 


tan 


Wn  VT^  DR)2 

'  c?-<DR,w* 


(1-29) 


wn  VT-IdrT2 

V»  ■  T7{ 


Practical  values  for  (DR),  W  ,  and  (CT)  are  given  in  part  B,  Eq.  (1-26).  If  these  are  substituted  in 
Eq.  (1—29),  the  result  gives  the  transient  response: 

[(C)V](t  i)  =  ( [(C) V](tfi))0  {1.078  e“  1-950  rio  4  t+0>18  e-354  x  10  4  »sin(.ooil26t 
-  0.437)  }  ♦  U(C)V](t  i))o  {  549  e  - 1-950  x  10  4  t 

♦  976  e  "  154  x  1°  ’*  f  sin  (.001126  t  -  a 594)  > 

♦  ([(C)V](  i))0  {  775000  e“1950  *  10'’  f 

♦  78100  e  ~  3-54  x  10  4  t  sin(.001126t  -  1.71)  }  (]_30) 
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D.  EFFECT  OF  j  AND  y  COUPLING. 


In  part 


A  if  it  is  assumed  that  i(jnt) 
^(IV)i^x 

*(  int)  x  e 

nfcmds)(i,V)x 

i 

^(cmdsXi.W)* 


~  '(int)x  aml  (^(gu)  = 

■5  '  .  -pMC)V)(1-i„) 

n(cmdsXi,^  )x 

{  [W,F  ♦  (Lim)t!  cos (Lot)t  -  p  [(C)Vj(t>i)x  } 


(1-31) 


This  applies  to  a  single-axis  system.  In  a  two-axis  system,  the  Earth  rate  couples  the  vertical  indicators 
as  follows.  Define 


'( int)  x 


(1-32) 


Then 


s, 


(cmds)(i,U’)  x 


{ [WIE  *  (Lon)t]cos(Lof)t  -  [WJE  ♦  (Un)t3  [(C)  V] (t  i)ysin  (Lot), 


-  P  KC>V)(tiI)s 


(1-33) 


where 


(Lon)t  cos(Lot)t  -  (Lon)t  [(C)V](t  i)y  sin(Lat)t  =  W(EV)t  * 


(1-34) 


An  analogous  situation  exists  for  the  y-system 

1 


-,(cmdsXi,«')y 


-<  -M  *  '"it  •(LW!,[(CIV](l  illSln(U),  -p[(C)V](li0,:i 


U-V) 


where 


-(Lot),  ♦  (Lon)t  [(OV](t  i)x  sin(Lrt)t  =  ^(FV)ty 


(1-36) 


The  x- and  y- systems  are  thus  coupled  through  two  terms,  each  proportional  to 
[WIE  +  (Lon)t]  sin  (Lot^.  for  purposes  of  the  coupling  effect  analysis  given  in  the  text, 

assume  that 


(Lon)  =  0] 

‘  >  i.e.,  the  system  base  is  stationary  on  the  Earth. 

(1-37) 

(Lilt);  =  oj 

(1-38) 

^(cmds)(i,W)x  ^(cmds)(i,W)y 

(1-39) 
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Insetting  tnese  values  of  it  and  iy  into  the  development  in  pan  A,  the  two  performance  equations  in 
z  and  y  (ignoring  uncertainties  in  component  outputs  as  well  as  inaccuracies  in  system  inputs)  are 

(p  *  fcV)  )(pJ  *  2(DR)WnP  *  "n>!(c>V3(t,i)y  -  { p 2  -  l2(DR)Wn 

♦  (cY)]p>wif  sin5-af)i  ^C)V](t>i)x  a-40) 

(P  *  2(DR)Wnp  eWB*)«C)V](tii)l  - 

-  {PJ  ♦  [2(DR)Wn  ♦  (C!r)  J  p  }  W1E  sintLofJj  [(C) V ] (ti|)y  U-41) 

It  is  on  these  two  equations  that  Fif.  7  in  the  text  is  based. 
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DERIVATION  2.  SINGLE-AXIS  POSITION  INDICATORS 


A.  OPEN-CHAIN'  INTEGRATION  OE  THE  ANGULAR  VELOCITY  OF  THE  INDICATED  VFRTICAL 


Refer  ro  Figure  2—1.  The  indicated  position  is  given  by: 

'(ixdo  4 


P(i)  *  S(P«XiP)  p  '(hfi)  *  p 


From  Eq.  (1-8)  in  Her.  I  A, 

1 

(cmdsXiJI) 

^(EV)i  *  ^(F.V)t  “  P  C(C)  V]  (t.i) 


Vint)  *5^1  tV/(Kv)i-(U)W(gu)-(l)W(,EXcp)] 


(2-1) 


(2-2) 


[1-9] 


Vi)(i,P)  *  S(cmds)(i,W) 

(This  constitutes  an  arbitrary  adjustment  of  system  parameters.) 


(2-3) 


Note  that 


and 


7*W  -p.-ip). 


-r>  ■  KC)VI(tJ)  - hov:,.,,,. 

From  the  foregoing,  Eq.  (2-1)  mey  be  expressed  as  follows; 

P,  •'  P,  -(Pt)„  -  [<C)V)(tii)  .  I(CIV)(,.ite 

-7«u)»(rl|  ‘(P).  *7(uiP(ixd0 

Define  the  correction  to  the  indicated  position  as  follows: 

«C)P](t(i)  *Pt-pi 

At  the  start  of  the  navigation  problem,  this  becomes 

«C)Pl(t>i)o  •  (Pt)0-  (P;)0 


(2-4; 


(2-5) 


(2-6) 


(2-7) 


(2-8) 
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Substitute  Fq.  (2-7)  and  (2-8)  in  (2-6).  The  result  is 

mcip:m  ■  t<CH*iw„  •  ttovjw,  -[<c)V)(lllo 


(2-9) 


B.  OPEN-CHAIN  INTEGRATION'  OF  THE  ANGCLAR  ACCELERATION  OF  THE  INDICATED 
VERTICAL 

Refer  to  Fig.  2—2.  The  vertical— indicating  function  of  this  system  has  been  analyzed  in  Derivation  1. 
The  position-indicating  function  is  analyzed  as  follows.  From  Fig.  2-2,  fas  in  Derivation  1,  Fq.  (1-3)  ) 


1  , 


1  ,...r 


f/i.r:v  =  S/Lt:w*  is  e( fl»)  5  *  ^J'(hfi)(dO  *  '(hfi)o 


*( hfi )  =s(hfi)(e,i)  p  le(dir>  *e(lfi)  ’e(cp)' 


'(hfi)o  *  5 


1 


'(cmds)(i,T) 


"(EV)s 


(2-10) 

(2-11) 


»here  ^jry)s  is  the  initial  setting,  in  the  second  integrator  in  the  vertical  indicator,  of  the  angular 
velocity  of  the  indicated  vertical  with  respect  to  the  earth.  Also,  from  the  figure, 


therefore 


e(fb)  ’  S(fb)(i,e)'(hfi)  *  (UKfc) 


-'l  *S(hfi)(e,i)  ^fbXi.e)  ~ ''(hfi)  *  ^hfiXej)  ~  [e(dir)  *  e(Ifi) 


♦  —  (U)i,.  .o  * 


e(cp)  ~'u,e(fb)J  *  -^u>'(hfi)(d0  *  C  W(EV)s 

F  (cmdsXi.W) 


a- 12) 


( 2-1?) 


Also, from  Fig.  3—2, 


*(dir)  *  S(dirXe,e)e(au)  +^U^e(dir) 


1  1 

e(!fi)  ’  J(lfiXe.e)  ~  e(au)  *  ~ 


*  Wv*  .d  —  e/...t  ♦  —  (U)e{!fiXdf)  +  *(tfi)o 


f2-M) 

(2-!<) 


The  constant  ^{fj)0  can  be  evaluated  from  Fq.  (2-10),  with  initial  conditions  inserted  such  that  the 
input  to  hfi  (high  frequency  integrator)  is  nulled: 

S(fb)(i.e) 

#(Ifi)o  *  - - 

^cmdsXi.T) 


^(F.V)s  ^(dirXe,e)e(au)o  ^*(dir)o  *(cp)o  *^^e(fb)o  ^  ^ 
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•substituting  bq.  (2-14)  and  (2-1^)  into  (2-H)  there  results 


Define 


1  *  5(hfi)(e,iAfb>(i,e)  ~Hhfi)  =  VfiHe.i)  ~  !S(d.rXr,e) 

*  ^lfiXe.e)  j’e(au)  ’  S(hfiKe,i>  ~  '»fcp)  "  ^Kfb) j 

*  S(hfiXe,i)  “  ^U)e(dir)  4  "!U)*(lfiXdf) 


W»\ 


♦cnidsXi,*) 


(FV)s  ^difXo.e)  ®(auk)  ^^*(dir)n  *(cp)o 


-  (U). 


4  ^{fbVo  '  *  —  (U)*(hf;Xdf)  4  e  ~~  V:v>s 
P  i(cmdsX»,«’) 


1 

(CT), 


S(hfi)(e,i)  S(fbXi,e) 


(U) 

1  ^IfiXe.e) 


,CT)‘ld>  W-.« 

The  gyro  input  is(cf.Fq.  (1-2)  and  ( 1—3)) 

'(Hfi)  =  7“  {W(EV)t~  P^C)V](tji)  -  (U)W(gu)  -  (l)W(1EKcp)  > 

5(emdsXi.T) 

Solving  for  *(au) .  Eq.  (2-17)  becomes 

VfiXe.i)  ^dirXf.e)  4  (CT)(,d)  ^  *(au) 

^,>W(lE)(cp)  3“S(hfiXe,i)  J  '•(cp)"*(cP)o_^)e(fb) 


4^)»(fb)0  4  ®(dir)  ^*(dir)o  S(dirXe,e)  "(auk- 
♦  — (U)*/if:v  jrv  }  ~  ::  —  —  W/t.,,4 


7'u'*(lfiXdf)  '  (rf)  5  7  (EV)s 

P  ltl,(lg)  5(cmdsXi,W)  P 

1  1 

-  -(UJUfi 


P  (hfiKd0  i  aw.,.  <FV>S 

P  T(cmds)(i,W) 


(2-P) 


(2-18) 


(2-19) 


(2-20) 


(2-21) 
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1  VdiXe.P)  1 

(UKlf.Xdt>‘"  "  “  (U)' 

^(hfiXe.i)  ^(dir)(e,r) 


fhli/df) 


12-20 


It  will  be  assumed  that  the  system  be  initially  aliened  under  conditions  of  zero  acceleration;  then 


"(au) 


o  '  S(auXa,e)9  ^C>  V](r.i>o  *  (U>e(au)o 


(2-21) 


1 


1 


^(cmdsXt  V)  ^(hfiXe.i)  ®(cp)  ^(CT)(|  (CT)(!^  ^  I  v'Xcp) 


(2-28) 


The  position  indicator  is  calibrated  by  setting 


VdiXe.P)  1  S(hfiXe,i)  S(dirXe,e)  S(cmds)(i.«) 


Substituting  the  above  into  Kq.  (2-2(5),  there  results 


(P 


1 


1 


KTI(IJ))(Pi  ■  Pt)  '  KT)(ld):(P'’»  "(P-!»  :  ■  ■  (CT)„k,  "  <CT)„d|  !|P'  "(p>)”1 

■W  -(c^,W*^-(cTjr/i>.- 
■'(Cf),,:, '  (c¥w*p  ["<E™'r>  (r  •<CT&'«CWW) 


^(CT)  +  S(Pd,Xe,P)  S(auXa,e)9  p  ^C>V^(t,i)o 


(1 


'fig) 

1  1 


(CT) 


H(U)W,  .  -(I)W, 


(Ig) 


(gu)  T  w"(IE)(cpl  ' 


_S(hfi)(e,i)  Sfcmds)(i,«')  p2  P  ^U^*(dir)  ^U)e(fb) 


*VdiXe,P)  p  ^UKau)o  S(hfi)(c,i)  ^cmdsXi.T)  ^  tU)*(lfi)(df) 

"  VmdsXi.'F)  p  ^U>'(hfi)(df)  *  (P  *  (CT)  * (  p2  W)P(di){dO 


*  '  W^diwrd 


(2-50) 
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Define 


■  p,  -pi 

(p  '  iCT)(IJ),[IC>P1"’‘)  '  fCT^1(C1P1<*'i»  ''(CT),,;,  ‘(Cf),,;,1  P  ,*<EV» 

** W(KV)(cp)^  *  (CT),,  )^W(KV)s*  “W(EV)(cp)o’) 

“  (CT)^((Pi)*#  +  (p  *  (CT)^  [{C)V]<m) 


(2-31) 


'(U) 


^(CT^J  J  *  ^(auXa,e)^ChfiX<v)^(dir>(e,e)S(cmds)(i,«,)9  p^  l(C ) ^] (tf|)0 


+  ^(hfiKe,i)^(cmds)(i,W)  t(U)e(dir)“  (U)e(dir)<>- (U)*(fb) 

+  ^UKft>)o  +  f  1  ^)i(hfiKdO"S(PdiXe,P)  ~(U)*(au)o 
*hfi)(e,i)  P 


+  ^hfiXej)  ^cmdsXi.W)  p2  >*(lfi)(df) 

-|p*,c^1)I^(u»W*7,U)p™1 


(2-32) 


Since  W(EV)s  and  (Pj)s  both  represent  the  best  information  on  the  initial  velocity  of  the  base  at  zero 
time. 


W(EV)s  *  (Pi>. 

Define  the  inaccuracy  in  the  initial  position  rate  setting  as  follows  : 

^W(EV)s  *  W(EV)s  "  ^(F.V)t^o 
and  the  inaccuracy  in  the  velov  y  compensation  as 

(WfEVX-.p)  3  W(EVXcp)  "  W(EV)t 


(2-33) 


(2-34) 


(2-33) 


From  Eq.  ( 1—9)  and  ( 1-13),  note  that 


VuXa.e^fhfiXe.b^dirXe.eAcmdsXi,*) 


(2-36) 


70 


(P  4  <CT)(k,))  [(C)P3(‘.i)  '  (CT),1U  ((C,P]«.')o  “(Ct),^  UC)V](t.i)o 


'(Id) 


'(lg) 


(p^CT^)[(C)V1(M>+{U)W(gU)*(«)W(M.Xcp) 

'  (ig) 


"(U)PW^n‘  l  {{(CT)(|f!)  "(CT)n,n)[(l)WC  v)(cp) 


'(id) 


+  ^W(in)(cp)'  *  ^hfi)(ff,i)  ^(cirds)(i,W)  '^U>e(dir) 


^^(dir)o"  ^)#(fb)  ‘(U)e(fo)o]  4  VmdsXi.VlJ^'fhfiXdf) 

p  c  ^^®(au)a  fJ)P(di)(aO 

s  lauXa.e) 

”(CTLTx  ^U)P(di)(dO  3  *  ~2  {S(hfiX<r,i)S(cmds)(i,«){U)*(lfi)(dO 
'  #(ld)  “ 

"fCT)(1J)  (U)P(^Xdf) 
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C.  APPLICATION  OK  DOUBLY-INTEGRATING  ACCKl.EROME'  R  TO  POSITION  INDICATION. 


Refer  to  Fig.  2—3.  The  vertical-indicating  function  of  this  system  can  be  analyzed  by  a  method  similar  to 
that  of  Derivation  1  to  arrive  at  a  differential  performance  equation  ill  [(C)V]^t  when  Schuler  tuning 
is  used  to  minimize  the  |erk  etror.  The  position-indicating  function  of  this  system  is  analyzed  as  follows. 
From  Fig.  2-3, 


'(dir)2 

S(dir)2(e,i)e(dir)l 

4  (U>(dir)2 

(2— 3ft) 

'(diff) 

=  S(diff)(e,i)Pe(dir)l 

4  ^)'(diff) 

(2-39) 

®(int) 

_  <  1 

5(int)(e,e)  e(dir)  1 

4  “^U)*(intXdf)  4  ®(int)o 

(2-40) 

#(dir)l 

^(dir)l(e,e)  *-®(dia) 

“  ®(int) '  4  ^)®(dir)I 

Define 


■*(dir)l(e,e)  ®(dia)  -*(dir)l(e,e)  ^(intXe.e)  .^®(dir)l 
“  S(dir)l(e,e)  —  CW)e(;ntXdf)  “  S(dir)l(e,e)  ®(int)o  *  ^U^e(dir)2 


*(ep)  "  {pF)(cpK«’,i)  W(EVKcp) 
e(dia)  "  ^(intXe.e)  ^ 


_ L_ _ 1, 

S  S  .  +  n  "  ®(dif)  1 

^(dir)l(e,e)  *(intXe,e)  p 


4  (U)®(intXdf)  4  e(int)o  “  c  ~  ^^e(dir)l 
P  5(dir)l(e,e) 


(2-41) 

(2-42) 


(2-43) 


(CT), 


(lg)  =  s  s 

5(dir)l(e,e)  3(int)(e,e) 


(2-441 


(CT) 


(id) 


S(diffXe,i) 

^(dir)2(e,i) 


(2-4*5) 


Operating  on  Eq.  (2-43)  with  (1  ♦  (C T )( id)  P-  an^  substituting  Fq.  (2-44)  and  (2—45),  the  result  is 


[1  ♦  (CT)(ld)ple(dia)  =  ^(intXe.c)  ^T) 


1 


1  (; 


1 


<»*>  P  S ..  1  "{dir)2 
K  J(dir)2ie,i) 


'(diff)  ^)'(dir)2  ^'(diff)**  *  S(intXe,e)^CT)(ld)fe(dir)l^ 


♦  [1  +  (CT)(M)p]  ^  (^)*(intXdf)  4  ‘1  4  (C1")(ld)P^  *®(int)o 

— ~ 

idir)I  (e,e) 


(2-46) 
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brom  Fq.  (2-41) 


^•(dir)l"o  “  ^)*(dir)lU  *  e(int)o 


®(dia)o  c 

a(dir)l  (e,e) 

Multiplying  Fq.  (2-46)  by  — ;z|r — *  and  substituting  1  0.  n-47),  there  results 

'CT'(ld) 

.  1  .  1  ,1  J. 

(P  +  (CT)(Id),,Wi®>  “  (CT>, +  (CT)fld)  (CT)f 


^(dir)l(e,e)  ViffXe.i) 


(ig> 

(1 


'(Id) 

1  in- 
(CT){I  ,  pJ  i,(dir>2 


'rig) 


'(Jiff) 


<u>'(dir)2  “  ^U>'(dirf)-  *{p  +  (cf)firip(U>*f,n,Kdf)) 


f  2-4?) 


S(dit)I(e,0 


(U)®(dir>I  +  S(int)(e,e)  KUKdir)l'o 


(2-48) 


At  t  ■  0, 


1 


W  +  L(CTy  -  -  {CT)  ]e(aia)o  -  l(CT)..d.  -  (CT)(1  ,  J'(im)o 


1 


1 


'(id) 


S(dir)l(e,e)  S(diff)(c,i) 


(Ig)  ’  '(Id)  "*  1  '(lg) 

{  t!(dir)2^o  *  '(diff)o  “  ^U)'(dir)2-o 


(U)i(diff)o>+  (U)®ri 


1 


(int)(df)o  V 

^(dir)l  (e,e) 


( (U)®fdir)l^o  +  *S(int)(e,e) 


(CTW  S(dir)l(f,c) 


)  ^U^e(dit)l^o 


(2-40) 


The  gyro  input  is 

'(gu)  =  i(dit)2  *  '(diff)  *  '(cp) 


'(gu)  £ 


'  -  -%:n  -  >  «c)»l„.i,  -  w)  v 


'(cmdsKi.W) 


(2-50) 

(2-51) 


At  zero  time,  it  is  desired  that 
W(KV)s 


and 


(g’J) 


(cp) 


3(cmds)(i,4) 


(2-42) 


(2-44) 


74 


75 


s 


'(indXc.P) 


[(1 


1  1 


Vdir)l(e,r)S(diffXe,i)S(cmdsXi,tt)  ^C1"W  P 


MW(FV)t  -P[(C)V!(tii) 


(U^(gu>  "  ^^W(lFXcp)  “  Vmds)'(cp)  “  W(l-V)s^  ~  £ 


S(ind)(e,P) 


*(dir)l  (e,c) 


{p(U)e 


!(dir)l 


-  ;(U)e(dir)l3o  ‘  (CD  3{U)e(dir)1}  -  .{U)etdjr)1;o  .> 


'(Id) 


1  1 

J(ind)(f,P)  l(CT)(ld)  p' 


m  ( >  (U)e(imXd0  ♦  (U)e(intXd^ 


(U)*(i 


(int)(df)o  ' 


S(ind)(e,P) 


1  1 


S(dir)Ke,f)  3(diffX^.i) 


i—  ^/CT\  -  ^U)'(dir)2  ♦  (U)'(diff)3  4  ^3'(dir)2 

5(diffX^.i)  '  '(lg)  P 


^U3'(dir)230  4  ^W)  _  (u3'(diff)o  3  (2~57) 

The  best  available  estimate  of  the  velocity  of  the  base  with  respect  to  the  Earth  is  Wfpv)s‘  Therefore 

(2-58) 


<P'i>o  =  W(EV)s  4  <Pi>(dt)o 

1  1 


1  .  1 


(CT)(1  p_)  S(cmdsXi,W)i(cP)  *  ((CT)(  "  (CT)(ld))p  W(^'Xcp)  (2-59) 


The  position  indicating  system  >s  calibrated  by  setting 

S(ind)(e,P)  *  S(dit)l(e,e)  S(diff)(e,i)  S(cmdsXi,V) 


then 


(p  +{Cf)(ld))Pi  =  (CT)MjPi)°  *  W<FV>S  4  (P^  4  W<F-VH 


1  1 


(CTV 


(lg)  P 


'(Id) 


’CKV)t  “(1  p) (p[(C)V](t  ;)  MU)W( 


1  1 

w 


(gu) 


^^(lEXcp))  -(CT) 


'(lg) 


(cirid^p  'Ev»r>  "<™* 


w, 


S(diff)(e,i)S(cmdsXi,'T)  3p(U3p(dir)l  "  ^U^*(dir)l  ^  o 


1  1 


(2-60) 


4  (CT)(H)  ^^*(dir)l  f(U)e(dir)l  3  p)  3  ♦  ^(indXe.P)  ^  (CT)f!4-j  p 


(Id)  P 


76 


4  *UfyirttXdf)  <U^(intXdf)o  *  "  VmdsXi.W)  *  (fT)  ^'<dir)2 

VV-  1  '(lg) 

4  {U)'(diff)’  4  <U^(dii)2  “  ^'(dir^o  4  ^U)'(diff)  ~  ^'(diftto  *  (2-fil) 


KOPl(tli)  -P,  "Pi 

il)  may  then  be  expressed 

1 


(2 — 62) 


in  terms  of  [{C)P](ti): 

(p  *  (CT)(ld))  [(C)P]('.^  =  (p  *  (CT\ldJ)  Pl  ~  (CT)(ld)(Pi)°  "  W(HV,S  ”(f5i)Mf)o 
-W<^-W<1  4(C^)7)(P[(C)V^> 

4  <u^w(gu)  4  (D^iEXcp))  4  ^(CT)(lg)  ~  (CT)^  P  W(EVXCP}  4  W(FV)s 

4  S(diffXe,i)  S(cmdsXi,«’)  *P^UKdir)l  "  ^U^(dir)l^o  4  (CT}^  ^UKdir)l 

t(U)®(cj£r) i ]  0  ]  }  -  S(indXe,P)  ^Trf)  ~<UKintXdf)  4  <UKintXdf) 

'  '(Id)  •* 

_  ^U)*(intXdf)o  }  4  S(cmdsXi,W)  {(CtT“71(U)W)2  *  ^U^'(diff) 

‘  "(1b)  P 


In  the  equation  above, 


7»(HV>.  '  P>  - 

(2-64) 

7W<EV*ep)  '  P,  "  |Pd<. 

(2-66) 

P  Pt  *  W(EV)t 

(2-66) 

Thus  Eq.  (2-6;)  n  ay  be  written: 

<>  •  left;-'  ■  '(cfi,;-  •  <ct)"~,(!<c>p1«!>>-  -  *w 
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•<>  -(CT),,,,  pHp!<C,V!".-)  *IU)V 

4  S(diHX^i)S(emdsXi,«)  (  p  (U)«(d,r)l  "  ^U)*(dir,l3o 

4  (CT)nd)  (U^(dif)1  ~  ^UKjk)»'o-  '  -  S(indKeJ')  (CT)(U)  p ^U,*(ini)<df) 

4  ^^(im)(df)  ^^(intXdf)o  4  ^(cmdsXi,*)  '  frT\  D  L^^*Cdir)2 

'  '(Ip)  p 

4  (U)i(dif0)  *  (U)i(dir)2  -  t(U)iwir)2J0  ♦  (U)i(dif0  -  (U)icdif0o>  (2-67, 
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DERIVATION  3.  AZIMUTH  STABILIZATION  SYSTEM 


Kefer  to  Fig.  3—1.  The  angular  velocity  of  indicated  north  with  respect  to  inertial  space,  about  the 
z-  asis,  is: 

^(l\’)i  ^(cmJs)(i,H  )z  'z  +  (U)W((t-„ 

This  is  related  to  the  angular  velocity  of  indicated  north  with  respect  to  the  Farth  as  follows: 

W(KN)i  *  W(!\’)i  ‘  V)sin(Lo,)t  ♦  tf(C)V)(t  i)yWir,cos(Laf)t 

From  these  equations, 

WUxiit).'.  ‘<U>"(gu)z  *WiE«-n(Ut)t  ‘  [(C)V]{t  i)ywIEc«riUt)t 

Define  i(C)N]  as  the  correction  to  indicated  north.  Then 

p iici'M] (t>j) . w(KN)t - w(EN)i - w(EB) .  iroHi,,., 


At  marine  speeds,  [Wpp  *  [(C)  N](,  *  0.  (See  Part  I,  pp,  BS-B6).  Therefore, 

p  {(C)  N]  tt,i)  *  W(KN)t  “  W(FN)i 

From  the  foregoing,  this  may  be  written, 

p  t(C)N]  (t  j)  *  W<EN)t -S(em<ls)(j,f)zi  *  *(U)W(su)z 

-WIE  sin(Lot)t  -  [(C)V](t  i)yW,F  cos(Lat)t 

From  Fig.  3—1. 

'z  *  {S(dirKi,i)z^y  *  '(cp)z>  *fU>'(dir)z 

“ix  (sin(Lol)j  *(U)Res)  Msec  (Lot);  +  (U)Res) 

where  (U)Res  is  the  uncertainty  in  the  resolver  generating  sin  (Lot);  and  S*C  (Lot);.  From 
geometrical  considerations,  i  and  ix  are  given  by: 


iy  *  r— — - -  {-(Lo»)t  ♦  [W,E  +(Lon)t]  [(C)V](t  ;)x  sin(Lot)( 

^(cmdsKi.V)y 

*  (W|E  ♦  (Lon)tl  ((C)N](f  i)co,(Lot)t-p[(C)V](tii)y-(U)W(?u)y  ) 


(3-1) 


( 3-2) 


(3-3) 


(3-4) 


(3-5) 


(3-6) 


(3-7) 


(3-8) 
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«et 


-  .  { lWIF  ♦  (Lon),]  eos(Lut), 

Tcmds)(i,tf)x 

-[WIF  *  (Lon),]  ((C)  V](,  J)y  sin(Lof), 


♦  (Lit),  [(C)H](,(i)*  p  E(C)V](,  j)x  -  (U)W(j|u)x  ) 


^cmds 


sXi,H)i  ^(cmds)(i,Vt')y  "  J(cmds)(i,»)z 


(3-9) 


(3-10} 


Combination  of  Eqs.  (3-7),  (3-8),  (3—9)  end  (3-10)  with  (3-6)gives  the  following: 

p  l(C)  N  ]  (t,i)  ‘  W(liN)t"S(dif)(i,i)z  '“(^t  +S(cmdsXi«)'(cp)zj  Ssec(Lo,)i  *  (U)Res. 

~S(dir)(i,i)Z[WlK  ♦<*•>,]  aC)V](,.)x(sin(Lo.),]  [ncMi  *(U)Re.) 
~S(dirXi,i)z[WlK  C(C) M ) (,ti) Eco*(Lot)t ]  [$ec(Lo.)j  ♦(U)Resj 

,Wi,i)t!sKH,(U)R«slP  [(C)  V]  (t,i)y 
*  S(dir)(i,i)z  [sec!Lo,)i  +  (U)Res](U)W(}tu)y 

+  [W,E  ♦  (Lon), 3  [cos(Lat),]  (sm(Ltjt);  ♦(U)Res]  [sec(Lat);  ♦(U)Resl 
- Wjp  sin(Laf)t  “(W,E  ♦(Lon),]  E(C)V]{,  i)y  [sir, (Lot),]  fsinJLof), 

+  (U)Res]  (sec(Lot)j  ♦(U)Res]  ♦  ((Lot),  [(C)N](t  i) -p [(C)V]{t  i), 

~(U)W(gu)x}  [ sin(Lat)j  ♦(U)Res]  (sec(Lot)j  +(U)Res] 

-  KC)  V]  (t,i)y  WIE  «s(Lot),  -  (U )  W(gu)7  -  S(CmdsX •  w)  (U)i(dif)2  (3-H) 


Define  the  latitude  rate  compensation  inaccuracy  as  follows: 

(l)(Lot)cp  ?  S(cmds)(i  T)i(cpjz-(Lot), 

Dt  ine  the  correction  to  indicated  latitude  as  follows: 

[(C)  Loti (,.!)  «■  (Lot),  -  (Lot) ■ 


(3-12) 


(3-13) 
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Then 


sinfLot);  *  sin(Lat),  -  l (C) Latl  ((1)  cos  (Lot), 
1 


sec 


(Lot); 


COS  (Let); 

s  sec  ( Lot)  t  “  ;(C)Laf](ti)  fan  (Lot),  sec(Lot)t 

Since  the  time  variations  in  (Lon),  and  (Lot),  are  small  when  compared  with  the  variation,  it. 
[{C)V](t,i)x,  [(C)V](ti)y  and  [(C)  N  ]  (ti) ,  it  will  be  assumed  that  (Lon),  and  (Lot),  are 
quasi-static.  This  simplifies  the  subsequent  analysis.  Note  that 

W(EN),  *  -  (Lon),  sin  (Lot),  -(Lot),  1(C)  V]  {t>i)x 
-(Lon),  cos  (Lot),  [(C)V](,i)y 

From  the  foregoing,  the  performance  equation  is,  ignoring  second-order  terms, 

tp  *  Wi,.>  [WIE  +  (Lon),]  -  (Lot),  tan  (Lot),]  t(C)N](ti> 

♦  {S(dir)(i,i)Zsec(Lat)tP  _tWIF.  +  (L°"M  ^«"(L<',)tsm(Lo,)t 

♦  cos(Lat),]}  [(C)V](t|i)y  -  ( tan(Lat),p  ♦  (Lb*), 

+  S(dir)(i,i)Z  [WIE  ♦  CLoo)t3 tan(Lo»)c >  [(C)  V](t  i)x 

>S(di0(i,i,2-c(Lat)t(U)W(gu)y-(U)W(gu)xton(M, 

“(U)W{Ru)z-S(dirKi,i)i(0(Lat)(cp)s«c(Lat)t 

-S(cmds)(i,T)tsec(Lo,)t](U)i(diD=t 

-  [WjF +(Lon),]  [sec(Lat),]  [(C)Lotj  (ti) 

♦  [(U)  Res](1  ♦  sin  (Lot),  cos  (Lot),)  [W,,..  ♦(Lin),] 

The  follow :ng  assumptions  are  warranted  if  the  max.mum  ship  speed  is  taken  as  40  Knots 

(Lon),  <<  WlF 
(Lot),  «  W,p 


Define 


S(dir)(i,i)Z  WIH  (CT)? 


( 4—14) 

f'-l*) 


r  3- 16) 


o-n 

(  ■  - 1-'1) 

<5-:o) 
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The  performance  equation  then  simplifies  to 

P  *  (CT)~  (CT), Wu.  lMclLol|,;P‘(C>V  (>. 

-W|F«e(Uf)t{[(C)V](tfi)v.  lOU  ^, 


i)y 


1  P  *  (CT)  [»an(Laf)t ]  :(C)V](t  i)x 


*(cf )? w„; 


(gu)v 


■,„(L..),(U)W(p0)J-(U)W(su|i 


“  Vmds)(i,V)  fsec(Lat)t]  (U)i(air)_, 

+  [1  ♦  sin(Lat)t  co$(Lct)t]  W1F  [(U)Resl 


M-21) 


Ignoring  x- coupling,  the  y- current  is 

'  «C) HI (t  0 W, R co*( Lot),  - p [<C > VI „  i>?.  qui*((u)y 


V  =.s, 


(3-22) 


(cmdsXi.W) 


By  a  procedure  similar  to  that  of  Derivation  1,  Part  A,  the  y  -  performance  equation  including 
z- coupling  can  be  shown  to  be 


<P3* 


1 


1 


(CT)(w/!‘rsp‘r,  (CT)(U)y  MIC,V;'..a. 

=  [p2  *  (CTT - p]  W‘H  I(C)N](t,i)C0s(Lo,)t 

1  '(lf.)y 

"  [(CT)(10y  ‘  (CT)n,n~  1  p(l)(La,)(cr) 


'(ld)y 


D2  *  ■ 


(CT), 


Hg)y 


p](U)W, 


<p  * 


(gu)y  DC  r  <CT) 

KsS(auXa,e)  1  >< 


)  (U)« 


(ld)y 


(au) 


"^(cmds)(i,W)  S(hfi)(e,'i>  ^U>efdirty  “^e(fb) 1 


"^(cmds)(i,W)  ^(hfiX",MUKlfi}(df) 
“S(cmds)(i,U)  P^U^'(hfi)(df) 


o-:3) 
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,  1  g  g  1 

Operate  on  Fq.  (  1-21)  with  [p5  ♦  p'  ♦  p  ♦  -  —  j,  utilizing  Kq.  (3-23)  to 

(CT'(lg)v  Rs  Rs  1  '(ld)y 

eliminate  ((C)  ^1  (t  ;)v-  The  result  is,  ignoring  *  -  coupling  except  for  x-gyro  unit  drift, 

(p4  +  (efi -  P'  *  ^  +  W2ik  ]p2 

(CT,(lg)y  Rs 


w  ^ 

■  '%  ‘  <CTI„d)y  ‘  SctT,  ’  *  (CTI,lg), lp  ■  Sf.fCt’^CTi;  ’  ,C)N'".'I 

■  tp’ .  (cfj^»J  *  g*  ’F,  (CT),,,;  [(C)Lo,l('.i> 

Sec(Lot>.  f  3  ,  _ 2^ _ 1 _ )  2 

(CT)ZWIK  P  (CT>ds)y  (CTW  P 

*  tjT  *  «!>,*  ,E  (rj.jy  “  -- )  •  P 

Ks  1  '(ld)y  1  '(lg)y 

•4-  7rr, - HI)(L«),  ,  -  tan  (Lot),  (p 5  *  Try, — p2 

R,  ICT)ndw  * 


^'i(cTw,lu,v,i*“c(U’,'i''' 


*  ‘(CT)^  *  R,(CT)/1E  lp  *  ^(CT)(ld,ytCT),W  ,E  1 
“(F'  *  (CT)(lgly p;  * i Tcf),,^’ |U)W<«^ 

sec  (Lot)  1 

— - (p!  -  t(CT)=W2  IF  .  r  - )p 

",W..riICT)»".E  ICT,"d)>' 

(CT)Z  u,2  S(cmdsXi,U')S(hfiXe,i)Sec(La,*t  f 

(CT)(ld)y  Ih  (CT)ZW1F 


-(CT)zW2ir.]p[(U)e(ditr(U)e(fb)l 


^(cmds)(i,W)  S(hfiXe,i)  “C^^t 

ctSe 


P “(CT)^  IF  HU)e( 
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3(cmisXi,«) 


sec 


(La»)t 


(CT)zWIh 


(p-(CT)zWIE)p(U)[(hflKd0 


•  (CT.-"1  'r/'imM,-' 


^)‘(dir)z 


’  0  *  sin  (Lot) ,  cos  (Lot),  )W  ,F  (p1  ♦  —  p‘  ’  p*  P  *  p”  ,rJ\ - )(U)Res 

(CT,(ll!)v  KS  KS(CT)( 


;(Id)y 


0-24) 


In  Derivation  1,  which  follows,  an  expression  for  (C)Lot  ((  in  terms  of  (C)V  ,(  |E  and  (C)N  ,ti| 

is  derived.  Solving  f!qs.  (i— 11),  O—  2D,  and  H— 21)  simultaneously  and  ignoring  -mall  terms,  thr.re 
results: 


pS-i-pt  JL  1  .  -1-)  D:  *  I  -  1_  l-l 

CCT)(,  ,  P  Rs  P  R.  V(CT)n,t  (CT) 


Vs  1  '(Id) 

1 


Rs  (CT),  '(CT)(ld) 

1 


+  (CT),  )p  *  Rs  ((CT)(ld)(CT),  (CT)Z>  (C)N  <«.»  “  '  i(CT),WIK  P‘ 


1  )  p3  + - 1 -  ( JL 


KT)zWie  '(CT)(lg)  (CT)(ld)  (CT)L  (CT)zV/j,.  Rs 


1  ->p’*tK  1 


(CT)L(CT)(l  ,  (CT)(ld)(CT),  Rs  (CT)2(CT)LWIh 


1  _  w  »  0  ,  9  1  (  1 

W.,  P  R. 


^T)(ld)^CT)Z"IF 


's  (CT)(,dy  (CT)(Id)(CT)zW„. 


-  WIE)j  sec  (Lot)  t  f(| )  Lot]  (cp)  -  ton(Lot)t  [p4  *  p'  *  P 


Ug) 

9  1 


i  icrw p’  w”w  •  *’  •  t  (ct>ov  r 


g  1  ,  1  1  g  1 

R;(CT)zW--(CT)(ld)  ^  (CT),  P  +  Rs  (CT)(!d)(CT),.(CT)zWn  (U)VW 


P4  4 


(CTW)P3  +  i  pJ+  fet \uf  {m«"' 
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sec  (Lot) 


«  ,  1  1 


0 3  ♦  { 

^  u i  r  v  / 


1 


R.W,  r  lCTIl'CT>,*lF  <CT)„s,<CT)2Wlr 

w, 


Wn°p‘‘  +  ((CT)nd'(CT)1  (CT),Wik  "(CT)(Ki))P  (U)efau’ 


"  S(crnd«Xi,w)  s(hfi)(e,i)  (C  T)  WIK  P  *  ((CT)[  (CT)ZW 


IK 


-  W1K)  pi  (U)e(lfi)(df)  -  S(cmdsXiiW)  sec(Lat),  L(ct)jrW1^"  P ’ 


(CT)l(CT)2Wie 


WIE)p3!  (U)i 


(hfi)fdf) 


^(cmds)(i,w)  ^(hfi)(e,i)  sec(Lat)t  '{CJ)^y  ^  +  ((CT),7CT)ZWIE 


-  WtK)p3]  [(U)e(dir)y  -  (U)e(fb)l  -  S(cm<js)(iiW)  sec(Lot),  tp*  +  p 


A  +  J-  1  - 
R,  P  Rs  (CT)(U) 


p]  (U)i(dir)z  +11  +  sin(Lot)t  cos(Lot)t]  WIK  [p‘ 


1  g  ,  g  1 


(CT), 


p’  +  p2 


dg) 


Rs  (CT), 


p]  (U)Res 


(■j-’M 


(Id) 


R6 


8? 


Then 


p  (C)L*’(til)  =  S(dir)(i  i)L  (I KLrt)((|l) 

*  wn,  !(C)N](t  l)  cos(Lot)t  -  P  1(C)V  (t>i)y 

-  }  -  (l)(U)fcp, 


Transposing  terms  in  Kq.  <W1)  (ignoring  t  -  coupling,  exci  pt  for  x -gyro  drift),  solving  (or 
[(ON]  in  terms  of  p[(C>N\tti).  and  multiplying  through  by  (CT)2Wjpy  gives  the  following: 

W,F  [(ON:(t  i)  =  -W.j.  tCT^pKOM]^;,  ♦  «(L«i)Ip[(C)V;.(t>i)y 
-<  (CT)zsec(La»)t  U(C)V((ti)y  •  i(C)U](t  i)} 

*  sec  (Lot),  (U)W(gu)y  -  W,h  (CT)z[tan(Lot)t](U)W(gu)r 

-  Wn,(CD,(U)»W  “  C)(Lo.)(cf)mc(Lc#). 

-  W1F(CT)Z  S(cmds)(i  LSec(Lat),  (U)'(dir)z 

♦  [1  +  sin(Lflt)tcos(Lat)1 1  Wjj  (CT)Z  [(U)Res  (4-9) 


Define 


S(dit)(i,i)L 

^(dir)(i,i)z 


WIH  '  S(dir)(i,i)lJCTh  WIF  ”  (CT)L 


Substitute  Fq.  (4—7)  in  F.q.  (4—6),  utilizing  Kq.  (4—8).  Ihere  results, 


tp  *  ,CT)L1  [(C,U*,e.‘> 


1 

(CT)l 


[(C)  V] 


(t,i)v 


■»1F(CT)L  p,,c,N1«-i>m(U1' 
-  ""Wop.  -  (CT1),  W,F 


^cmdsXi.Vf)  . 

(CT)l  W„; 


(4-10' 


♦  (1  *  sin(Lrt)t  cos(Ljt)t]  [(U)Res]^..^  cos(Lat), 


R8 


(CT),  W„  C0S(U,)-  (U)V» 


i  4-11 


ro  eliminate  (C)V  ^  anil  (C)N  ^  l)t  I\q.  M-21)  and  H— are  utilized,  ll»e if  simultaneous 
solution  is  (oiiiittiug  small  quantities!: 


ps  *  J.  p4  *  p’  *  ~  (-  1 


(CT) 


(Ijj) 


1 


K  K  (CT)(U1  (CT), 


)  v2 


9  (.  1 


1 


Rs  (CT)(1J)  (CT),  (CT),  (CT), 

9  1 


)  P 


'  R„  (CT){Id)(CT),  (CT)/  (C)Lo'  *  (CT),  (CT),(Wn)-  p< 


1  (-  1  — 1  _»  * 

(CT),  (CT) .(W,,.) i  l(CT)(lftt  (CT)(ld, 


g_  i  3  r  i 

Rs  (CT),  (CT)”  (W,K)J  P  *  (Ct)„M1(CT),  (CT), 


9  ,  1 


1 


Rs'(CT)n<li  (CT),  ^  ‘  rV(CT)(H  (CT), }  (l)(Lo,) 


(cpl 


1 


(p 


1 


Rs  (CT),  (CT),WU,  (CT)(IlJ> 

9  1  1 


)  sin(Lot);  (U)W(f,uk 


,if  ‘  (CU",(U,W'«“'> 


Rv  (CT),  (CT),(W„.)J 
t  <CT),('cT),»„  ,p  ■  jcfj^)»K(U)V„ 


_L_  ,_L_  1  _J _  , 

(ct),  (ct),  l(w,,r  p  +  (CT)f,d1(w„.r  p  * p 


■)rr^ - (!!)«(,„) 


(CT)(|d)  RsSfaaVa.r) 

^cmdsXi.Vn^hf.Xc.it^j^^j^^^y?  P’  ‘  P)l(U)e(dirH  ~(U)e(fh>] 
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’  S(cmJs#«,1DSfhti)(e,il  (CT),^ (CT),  ((W1K)5  P  ' 

*  VmdsXi.rn  (C7)t  (CT),  ((Wir)J  P  *  P^^'fhf.xdn 


1  :(  1 


1 


♦  V.mi,*..*,  {CT)i  wu  '(CT^  (CT),  >P  R,  (CT), 


g  1 


g  1 

Rs:  (CT)fid)  (CT), 


]  (U)i(Jir)7 


♦  -X.  — c&s(Lot)t  ^  sin(Lo»),  C0f.(Lat)t]  (p  ♦  .  )  J)Res  f4-12i 

R,  (CT),  (CT),  1  *  (CT)ad) 


In  Derivation  9  of  Hart  l  it  was  shown  that 

[(C)Lot](tii)  -  r(C)V](ti)y  -  UC)PA  (til) 

where  [(C)PAl(t  ,,  =  correction  to  the  indicated  polar  aius.  [(QVl(t  j)y  becomes  of  interest 
in  the  discussion  of  roll  and  may  also  be  obtained  through  the  simultaneous  solution  of 
Eq.  (3-21),  (3-23),  and  (4-11),  as  follows: 


(4-13? 


<P'  * 


1 


4  ?3  S  •  1 

a«  a.  — n J  A  - 


(, 


’  1 


(CT>„„"  R/  R,  (CT),ld,  (CT),'"  R,'(CT)„d,(CT), 


]  »  _9_  _ ) _ _ \  (C)V' 

(CT)l(CT)/)P  Rs  (CT)(ld>(CT)L(CT),  J(,'n> 


{( 


1 


1 


1 


!CT)„„  '  (CT),  (CT)„,,)p'  *(CT),  ((CT)7  '(CT>(IJ)!  p’ 


’  (,’ 


1 


•  (cTwETM^H,KLi"W'  •  v  •  <(CT K 

*  (CT)(|;|)P’  *  (CT)L(CT)„„  Pl  ‘  '(CT),  *  (CT)„fl)p’ 

*(CT)7((CT)'  •  Fn“,P'  '(CT)l„j(CT)l.(CT),p)(U)"(e": 


00 


Wl|.  cos(Lot)t[p!  *  ( 


i 


Rs^uXa,cJP  ’((CT)(U,  *  (CT),)p  *  (CT)7  ((CT),  ’  (CT)(lij;,p 

*  (CT)(KJ/CT),(CT)J(U)4'aul  ~  S<c.uN)(i,»>S(,lfl*(ef;tlpJ  4  (CT)*P 

*  (CT)|(CT)J(U>*n,iKdn  ‘  *  (CT)_  P 

(CT)[(CT)?  P^^^l(1,Hydf'  *  ^ cmd^ K 4  P* 

*  (CT),  (ClT,  P^^U^*fdit>y  “  ^U^effb>3  4  SfcmdsXi.w)W|l  V 


'(CT)0g)  ’  (CT),  )p  ’  (CT)flK)(CT);  P 

-  W,3f  cos(Ltrt),  [l  ♦  sin(Lo»)t  cos(Lot),]  [pJ  *  (—|~  *  — )p 

(LI).  (CT)|i.| 


P  (U)i, 


(CT>n^)(CT), 


p]  (U)Res 


!4-l4) 


01 


DERIVATION  5.  LONGITUDE  INDICATION 


Rrfer  to  f  ig.  5_j.  The  longitude  titive  (upper  right-hand  cornel  of  the  figure)  rotates  at  the  indicated 
celestial  longitude  rate: 


p(Lon). 


(cel)i  S<ldsKe,'J  )  *(sg) 


(5-1) 


where  (Lon)(ct|)-  is  the  indicated  celestial  longitude  and  is  the  output  of  the  integrating  gvro 

unit  mounted  on  the  longitude  gimbal. 


(L£n)(cel)i  =  W(IF>i  *  (Mj 
*(**)  *  *(sgXA,e)Ag 

where  is  the  gyro  gimbal  angle.  Define  Wjp^  as  the  angular  velocity  with  respect  to  inertial 
space  sensed  by  the  longitudinal  gyro.  Then  the  torque  balance  for  the  gyro  unit  is 


(5-2) 

(5-3) 


(5-4) 


(5-51 


(5-6) 


(3-1) 


HiW(PA).  ♦  (U)W(j?u)I  -  p(M(cd)i]  *  cpA^ 

From  the  foregoing 

P  ^(cel)i  ’  ^IdsKe.lD^sgXA.e)  [W(PA)i  *  <U)W(gu)l  ~  P(Lfln\cel),^ 

^(PA)i  mil^  ^  evaluated  as  follows.  From  geometrical  considerations 

W(PA)i  *  {W(iv)i)*  cosfLtt);  -  W(1N)i  sin(L*); 

From  Derivation  3,  Eq.  (3—1), 

WdN)i  ^(cmds)(i,T) *z  *  ^^(gu)z 

Substitution  of  F.q.  (3-16)  into  (3-6)  and  solving  for  S(cm{js)(i  T)  i  z  results  in  the  following 

W(iN)i  *  -  (Lon),  sirt(Lflf\  -  (Li)t  [(C ) V ] ~  (Lon),  feosduf),]  [(C) V]  ;t  i)y 
p  [(C)N]^t>jj  -  WjF  sin(Lot)r  -  [(C) V] ^jp  cos(Lo»)r 

^avnK  =  ^cmdsXi.T)'*  *  ((J)Wfguh 
Solving  F.q.  (5—0)  for  VmdsXi.W, *«  *  <U>W(gu)z-  there  results 

<W(fV)i>K  1  Wif  *  (Lin),]  cos(LtH)t  -  [W]K  ♦  0in)t]  [(C)V](f  ;)v  «n(L*)t 

[(C)N]((>i)  -p[(C)V](t  i)t  (5-0) 


(5--) 

(5-8) 
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From  Eqs.  (5-  7)  and  (5-9),  (5-6)  becomes 

W(PA);  ---  [Wj,r  ♦  (Lon)( )  IcosfLoOj  cosCLof);  ♦  sirvfLrt),  sim’Lot^  ] 

+  U-bt)t  [(C)N]({ co^Lcrt);  *  {(Ln^sinlLofy  -  [cosOjoO;]  pi  [(C)V](()1 
♦  [sin(La»),]r{(C)H](t>i) 


Define 

t(C)U](tii)  MLa»)t  -{ut\ 

Then  in  the  first  term  on  the  tight  in  Eq.  (5-10),  assuming  [(w)Laf]  ^  is  a  small  angle, 
cos(La»)t  cos(Lrt)j  ♦  sin(Lat)t  sin(Lat)j 

■  cos  [(Lof)t  -  (Lat)j3  =  cos  [(C)(Lat)]  (t  ^  3  1 

ignoring  second-order  terms,  Eq.  (5-10)  now  becomes 

W(PA)i  *  *IE  +  ^L°n^t  +  (L°,)t  [(C)N](t  ;)  cos(Lot); 

♦  RLit)t  sin(Lol)t - [cos(Lnt)t]  p]  [(C)V](ti)x 

♦  [sin(Lot)t]  p  ['")N](t  i) 


Eq,  (5—5)  may  now  be  written 


,H, 


[P  *  S(ldsXc,W)  S(sgXA,e)  ]  P(L<>n)(cel)i 

,  H 


=  S(lds)(e,U')S(sg)(A,e)(c) 


WIE  *  (Lon)t  ♦  (Lbt)t  [(C)N](r  ;)eos(Mt 
♦  {(Lot)t  sin(Lat)t  -  [cos(Lot)t  ]  p  }  [(C)V](t  i)l 
+  lsin(Lot)tj  p[(C)N](tji) |  *  S(lds)(c >T)S(3g)(a,e)'  c  ^U)W(gU)i 


Define 


(CT)e 


S(lds)(e,W)S(Sg)(A,e)(c') 


Note  that 


(Lo">(cel)t  "  W  IE 


♦  (Lon)t 


ii 

(5-10) 

(5-11) 


(5-12) 


(5-1?) 


(5-14) 

(5-15) 

(5-16) 
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fp  *  (CT), 1  p(Lfln)(«i)i  ■  (Cf,  -  p(Kcel),  4  ^)tcos(Ltt)t 

♦  [sin(Lrt)t]p  >  UC)N](t  i}  ♦  {(Lrt)t  sin(Lc»)t 

-  frWi*)  w  nnyi,  ♦  1  (u)w, 

-  'i  ■  -  It, us  (CT)  (gu)l 

Define  the  sidereal  time  drive  angular  velocity  W(tjs)  as  follows 
W(tds)  *  P(L"»)(cel)i  "  p(Lon)j  =  W(IF); 

Substitute  Eq.  (5—18)  into  Eq.  (5—17)  and  multiply  by  (CT), 

[1  ♦  (CT ), p]  pd-on),  ♦  (1  ♦  (CT),P)W(tds)  .  p(Lon)t  ♦  W,E  ♦  (U)W(gu)1 

-  [(U)t  cos(Mt  ♦  (sin(Lrt)t)  p]  [(C)H](tfi)  ♦  [(Lbt),  sin(La»)t 

-  (cos(La»), )  p]  [(C)V](t  i)x 

Define  the  inaccuracy  in  the  time  drive  angular  velocity 
0W(tds)  «  W(tds)  -  WIE 

Define  the  correction  to  indicated  longitude: 

[(C)Lon](t  i)  =  (Lon)t  -  (Lon), 

then 

(1  ♦  (CT), p]  p[(C)Lon](t(i)  -  (CT),pJ(Lon)t  ♦  (1  ♦  (CT),  p)(l)W(t<!#) 
♦(CT),Wie  -  (U)W(fu)1  ((Lot), cos(Lrt)t  ♦  (sin(La»)t)p]  [(C)N](ti) 

-  [(Lif)t  sin(Lat)t  -  (cos(Lof)t)  p]  ((C)V](tj)x 

The  variation  in  Earth  rate  is  very  small;  therefore, 

pW,E  *  0 

Integrating  F.q.  (5-22),  there  results 

[1  *  (CT),  p ]  [(C)Lon](t  i)  -  [(C)Lonl(t  ;)o  *  (CT),  [(C)LW(t >i)o 


(5-17) 


(5-18) 


(5-19) 


(5-20) 


(5-21) 


(5-22) 


(5-2?) 
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*  (CT),  [ p(Lon)t  -  (p(Lon)t)0  3  ♦  (CT),  UDW(fds)o) 

*  P  " y  cos(U)t  ♦  (sin (Lot),)  p]  i(C)N](t  i) 

-  [(La»),  sin(Laf)t  -  (eos(Lo»)t )  p  s  l (C)  V  (t  l)x -f  (5-24 
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DERIVATION  6.  GROUND-SPEED  INDICATION 


1'he  true  ground  speed  is  i  where  is  a  set  value  of  the  Tarth-radius  and  W(  j  is 

given  in  terms  of  irs  components  hy: 


W(FV,t  ’  V  (Lon)t  cqs(Lq*)(  2  +  (Lot),'  U-'i 

From  the  above,  the  true  ground  speed,  which  is  the  magnitude  of  the  true  velocity  of  the  base  with  respect 
to  the  Karth,  is: 


r(HH}t  “  Rs  V(Lon)t  cos  (Lot),]  2  +  (Lot),2 


<(-2, 


Oase  I.  The  vertical  indicator  as  a  source  of  ground-speed  data. 

'he  quantities  in  Kq.  ( fs — 2 )  are,  in  part,  implicit  in  the  x-and  y-gyro  units*  orientational  control 
signals,  as  follows: 


W(IVXi)x 

^cmds)(i,«) 


- t  f(W,K  +  (Lon),:  cos  (Lot), 

i(cmdsXi,«) 


W,K  *  (Lon), j  ;sin (Lot),]  [(C) ♦  (Lot),  f(C)N] (t>;) 


~  p r(C) v j (ti)x  (u)w( 


"(gu)x 


(6-v-, 


w, 


nvxijy 


^(cmdsXi,*')  ^(cmdsXi,'*) 


-  (Lot), 


♦  [W[{r  ♦  (Lon),]  [sin(Lot),]  [(C)V](,  i)x 

♦  fW,R  ♦  (Lon),]  [cos (Lot),]  [(C)N](,  -  pi(C)V](tA  -  (U)W(gu)>)  <A-4> 


Define  the  gyro-unit  signal  compensating  for  the  horizontal  projection  of  the  Farth ’s  daily  rotation  with 
respect  to  inertial  space  as  follows: 

wapxcp) 

'«•  s  ~ 

J(cmds)(i,W) 


cos  (Lct)| 


o'— ^ 


The  quantities  lv  and  may  be  combined  as  follows  to  obtain  the  indicated  ground  speed,  in  a 
way  analogous  to  Kq.  (6-2): 


(KV)i  *  VmdsXi,*) 


in) 


*  t. 


(C—n 
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From  1-q.  (6-1),  (6-4)  and  (6-6), 


r(EB)i  -  Rs  •  (Lon),  cos(Lot),  •  W,,.  cos  (Lot),  -  W(])  )fcp }  cos(Lof)1 

-  W[f  *  (Lon),  sin(Loi),  (C ) V  f , (1  >y  '  (Lot),  (C)N  i 

-  P  (C)v’(t  i)5t  -<u)w(£U),-  *  ■  *  w!« 

•  (Lon),  sin(Lot),  ’(C)V  (,;)x  •  W|(  *  (Lon),  cos(Lot).  (ON  (t  ;, 

-  p-(C)v:(t(i)>  -<u)w,buI).]’>’!  (r” 


Define 


<I>WP)  W(IFKcP)-WU- 


((S-fel 


and 


(C)Lot](t>i)  -  (Lot),  -  (Lot), 


From  Kq.  (6-0), 

cos  (Lot)*  *  cos  (Lot),  *  (C)Lot  sin(Lot), 

Substitute  the  last  three  relationships  into  Hq.  (6-”),  and  expand  the  square  root  by  the  hinomjna. 
theorem,  utilizing  Kq.  (6—2)  for  ^0te  also  th. 

R  (Lon)  cos(Lat) 

— -  - -  cos(Az), 

v(F!H)t 


and 

-  R  (Lot) 

- — - 1  -  sin  (Az), 

v(F.H)t 


The  expansion  of  the  square  root,  ignoring  second-order  quantities,  then  results  in  the  following: 

v(t:n)i  '  v(hbh  *  Rs  ^  ”  (|)WdHXcp)  "  wn;  ^c)Lot  (t.n  sm(Lo,)t 

-  wn.  +  (Lon),'  sin  (Lot),  (C)V'(,iH.  ♦  (Lot),  ;(C)N ;  (t  j 

-  p[(C)V](ti)j,  ~  (U)WfK,0x)cos(Az),  *  (  W,,.  *(Lon),]  sin(Lot),  [(C)V-(l 


,i)x 
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♦  Wlt  •  (Lon),  cos ( Lai) t  (C)N  (t  lJ  -  p  ( C ) V  f[  i,, 

-  (U)  W,  ,u  v  )  sin  (Az), 


(  .»se  II.  '^elt-erccting  latitude  indicator  and  ce!esiia!-F 'art  h*rate-tr«t<.  King  longitude  indicator  as  sourer s  of 
ground-speed  data. 

l’hesr  indicators  are  discussed  in  Derivations  4  and  \  rcspectiveh.  My  differentiation  of  I  q.  (4—0  in 
Derivation  i,  the  indicated  I  at  i  rude  rare  can  be  given  by: 

pfLot),  -  (Lot),  -  p  (C)Lot  ([i|)  -  Sn  ux;t.r  i, 

From  equations  {*— M>)  and  (*—20)  of  Derivation  €t  if  (CT)j  is  sufficiently  small, 

pfLon),-  "  (Lon).  -  (l)WtN  *  (U)W(i:u)|  *  (Lit),  ; (ON  ft  cos(Lot)t 
♦  sin(Lat)t  p  (C ) N ; (t  i)  ♦  (Lot),  (C)V.;{t  i)x  sin(Lat), 

-  cos  (Lot),  p  (C)v;(ti)x  -  s(1;jsXl)#)  i,  <*-U) 

^ or  ground  speed  i n J i c ar ion  using  the  foregoing,  consider  thF  following  expression  for  vr p  p  1  i 

vn-:mi  ’  Rs  vU>(Lonjt]  cos(Lot);  *  p(Lat)|’ 3  (6-:r-) 

These  quantities  are  obtained  from  the  system  as  follows.  Longitude  rate  is  obtained  as  a  current  derived 
from  the  output  of  a  tachometer  geared  to  the  longitude  indicator.  If  the  sensitivity  of  the  tachometer  and 
its  associated  circuits  i<  Sy,j  and  the  output  is  iy, 

p(Lon),  -  .  ~  i1  (('-}') 

*rn,i) 


Indicated  latitude  is  proportional  to  tor  latitude  gimhal  drive  signal,  from  Fq.  (4—2)  of  Derivation  4: 


P(L°»>i  ’  ~  Sdds)  (i.vr)  h. 


(4-2) 


From  the  foregoing, 

Vf  F.B)i  *  ^ 


iT  cas(Lat);]  3 


*  S(lds)(i,U)  *!.. 


(6- IS) 


The  square  root  may  be  evaluated  as  follows.  From  Fo.  (6—10): 

cas(Lat):  ^  cos(Lat),  *  (OLtrt  ^,  ^  sin(Lat), 


(6-  in) 
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so  that,  from  eq.  (6—1*)! 

ip  (Lon)  cos  (Lot)  -  (Lon),  cas(Lot),  ♦  (Lon),  (C ) Lot  .(!  j)  sin(Lat), 

(I)  W  ( js  cas(Lot),  *  (U'W{^u.j  cos^Lot),  *  (Lot),  (ON  (ti ,  cos  ‘(Lot). 

•  'sin(Lot),  cas(Lat),:  p'(C)N.  ((1)  •  (Lot),  ((C)  sin(Lot),  cos(Lot), 

■;r- -}<-!) 

^  rom  Hq, 

v{FB.,  -  Rs  (Lon),  cas(Lot)t ;  *  *  (Lot),*  f'  -:i 

Therefore,  i*tnorin£  small  quantities  in  the  binominal  expansion  of  the  jqu.tre  font, 

V(KH)i  *  VU-  lilt 

♦  Rs  •  (Lon),  (C)Laf  (t<j)  sin(Lot),  -  (I )  W, ,  cos(Lat), 

♦  (U)W(tttt)1  cos(Lot),  *  (Lot),  (C)N  f,  ;,  cas2(Lot), 

*  ( sin(Lat),  cos(Lat),  p  ;'(C)  N ! (t 

*  (Lot),  '(C)  V ]  f t  i)x  sin(Lat),  cos(Lot), ]  cos  (Az), 

-  ;pf(C)Lat]f,ti)’  sin  (Az),  >  (  -20) 
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